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We  report  a  systematic  study  of  the  first  and  second  hyperpolarizabiliiies  of  several  small 
molecules  at  a  consistent  level  of  tlieory  and  basis  sets.  Coupled  cluster  (CC)  methods  forj 
correlation,  analytical  high-order  time  dependent  Hartree-Fock  (TDHF)  theory  for  I 

dispersion  effects,  and  polarizability-consistent  basis  sets  give  agreement  to  about  10%  between! 
the  calculated  hyperpolarizabilities  and  the  gas  phase  experiments  for  the  nonpolar  molecules,  f 

Hi,  Ni,  CO2,  and  C2H4,  and  effectively  nonpolar  CO.  Results  for  the  polar  molecules  FH.  i - 

H2O,  NHj,  and  H2S  are  improved  by  adding  lone-pair  basis  functions.  For  HjO  and  NH,  ^  - 

there  is  good  (  ~  10%)  agreement  with  experiment.  However,  a  --20%  difference  between 
experiment  and  theory  for  the  FH  molecule  persists;  this  difference  is  discussed  in  some 
detail. 
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INTRODUCTION 

Primarily  because  of  their  role  in  the  nonlinear  optical 
properties  of  molecules,  hyperpolarizabilities  have  recently 
received  much  attention.'’^  Hyperpolarizabilities  are  de¬ 
fined  as  the  coefficients  in  the  series  expansion  of  the  en¬ 
ergy  (or  dipole)  in  the  presence  of  an  oscillating  electric 
field.  In  the  static  case  this  expansion  becomes  £'(F)=£o 

-I-  •  ••  (Einstein  summation  convention  assumed),  with 
and  Yijiti  being  the  components  of  the  first  and  second  hy¬ 
perpolarizability  tensor.  In  addition  to  potential  applica¬ 
tions  to  optical  material  design,  which  requires  a  knowl¬ 
edge  of  the  P  and  y  tensors,  there  are  also  a  number  of 
other  unsolved  questions.  Some  of  these  concern  inconsis¬ 
tencies  in  the  experimental  values  which  cannot  be  ex¬ 
plained  by  the  different  dispersion  effects  encountered  in 
the  various  experimental  procedures.  However,  there  are 
also  significant  inconsistencies  between  theory  and  experi¬ 
ment,  and  between  theory  and  theory.  An  understanding  of 
molecular  nonlinear  optical  behavior  requires  resolving 
these  existing  differences.  The  objective  of  this  paper  is  to 
present  a  systematic  study  of  theoretical  calculations  for 
several  small  molecules  whose  gas  phase  experimental  val¬ 
ues  are  well  known.^’  This,  combined  with  comparisons 
with  other  calculations,  provides  some  calibration  of  the 
accuracy  of  high-level  quantum  chemical  methods. 

The  reliable  ab  initio  prediction  of  molecular  hyperpo¬ 
larizabilities  requires  several  elements:  basis  sets,  electron 
correlation,  frequency  dependent  (dispersion)  effects,  and 
vibrational  corrections.  First,  basis  sets  have  to  be  ade¬ 
quate.  The  prediction  of  hyperpolarizabilities  depends 
upon  products  of  matrix  elements  of  the  electron  position 
operator,  f.  Consequently,  unlike  the  molecule's  energy, 
which  primarily  depends  upon  inverse  powers  of  P,  any 
study  of  hyperpolarizabilities  must  allow  for  an  adequate 
description  of  the  more  diffuse  regions  of  the  molecule’s 
wave  function.  Basis  sets  that  are  usually  found  to  be  ad¬ 
equate  for  a  dipole  moment  require  further  extension  to 
adequately  account  for  the  polarizability  tensor,  a,  and 
even  further  extensions  for  the  first  hyperpolarizability,  P, 


or  for  the  second  hyperpolarizabiliiy,  y.  There  are  also 
other  factors.  Namely,  all  contributions  to  a  have  the  same 
sign,  so  the  larger  the  basis  generally  the  better,  and  this  is 
usually  true  of  y  as  well;  but  the  relative  values  of  the 
charge  density  must  also  be  described  in  p  (or  n)  and  that 
imposes  additional  restrictions  on  a  balanced  description. 
To  emphasize  the  basis  set  aspect,  the  SCF  values  for  static 
;'ii  =  +  (summation  convention)  for 

ethylene  are  —37  a.u.  at  the  STO-3G  level,  -241  for  a  DZ 
basis;  6230  for  a  6-3 IG-f  diffuse  (s,p)  basis,  and  6440  with 
a  (5s3p2d/}.s2p\d)  basis.  Of  course,  the  SCF  result  reflects 
just  a  part  of  a  basis  set’s  effect,  as  the  introduction  of 
electron  correlation  requires  a  much  larger  and  more  flex¬ 
ible  Hilbert  space. 

After  basis  sets,  the  accurate  inclusion  of  electron  cor¬ 
relation  is  the  next  most  important  factor  in  attaining  re¬ 
liable  values. The  initial  studies  of  correlated  molecular 
hyperpolarizabilities^’^  already  demonstrated  that  for  the 
small,  polar  molecules  FH  and  H2O,  correlation  could 
change  the  predicted  P  static  hyperpolarizabilities  by  as 
much  as  a  factor  of  2  and  have  a  large  effect  on  y.  This  fact 
has  now  been  supported  by  several  other  correlated  stud¬ 
ies,''’*  and  also  applies  to  recent  results  for  larger  mole¬ 
cules  like  para-nitroaniline.’'*  For  nonpolar  molecules  like 
/zanr-butadiene,”’’”  benzene,'^  and  the  nearly  non¬ 
polar  CO,'^  there  are  similarly  significant  (40% -60% )  in¬ 
creases  in  Pij  due  to  effects  of  electron  correlation.  Power¬ 
ful,  coupled-cluster  (CC)’”  and  many-body  perturbation 
theory  (MBPT)^'  methods  are  used  here  for  the  treatment 
of  electron  correlation  for  static  hyperpolarizabilities. 

Since  hyperpolarizabilities  require  third  and  fourth  de¬ 
rivatives  of  the  energy  with  respect  to  a  field,  a  purely 
finite-field  approach  requires  several  energy  calculations  to 
determine  all  tensor  components  and  can  sometimes  suffer 
from  numerical  imprecision.  By  virtue  of  computing  the 
induced  dipole  moment  in  the  presence  of  a  field  analyti¬ 
cally,  via  recently  developed  analytical  first  derivative 
methods  that  introduce  a  “relaxed  density, which  is 
now  i.-’p'c  ’TCiucd  f  r  ;hi_  CCSD  method  and  various  or¬ 
ders  of  MBPT^’  including  those  that  include  triple  excita¬ 
tions.  MBPTfTl.^"  CCSD  f  T(CCSD).  CCSD(T).  and 
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QCISD(T),^'  the  numerical  precision  can  be  improved  by 
a  couple  of  orders  of  magnitude. 

The  third  important  factor  is  the  frequency  depen¬ 
dence  0/  the  hyperpolanzability.  First,  without  the  inclu¬ 
sion  of  the  frequency  dependence,  the  various  experimen¬ 
tally  distinct  nonlinear  processes  that  occur  would  coincide 
in  the  static  limit.  For  /3,  i.e.,  second  harmonic  generation 
(SHG);  optical  rectification  (OR);  and  the  electro  optic 
Pockels  effect  (HOPE);  and  for  y;  i.e..  third  harmonic  gen¬ 
eration  (THG),  dc  induced  second  harmonic  generation 
(dc-SHG),  intensity  dependent  refractive  index  (IDRI), 
electro-optical  Kerr  effect  (EOKE),  and  dc  induced  opti¬ 
cal  rectification  (dc-OR).  Second,  dispersion  effects  can 
have  a  relartively  large  effect  on  the  observed  hyperpolariz¬ 
ability,  and  predictive  calculations  will  frequently  require 
more  accurate  dispersion  values  than  those  estimated.  This 
would  be  particularly  true  if  the  exciting  frequency  (<y) 
becomes  close  to  a  resonance,  and  for  SHG  a  resonance 
will  occur  at  (0/2  and  for  THG  &>/3,  further  compounding 
the  problem.  A  knowledge  of  the  frequency  dependence  is 
also  important  in  identifying  possible  resonance  enhanced 
effects  that  might  suit  a  particular  molecule  to  a  specific 
application. 

In  order  to  describe  the  frequency  dependent  effects, 
we  have  developed  an  open-ended,  analytical,  high-order 
time  dependent  Hartree-Fock  (TDHF)  approach^*  for  hy¬ 
perpolarizabilities.  This  permits  the  analytical  evaluation 
of  any  order  hyperpolarizability  tensor.  In  the  static 
((u=0)  limit  this  reduces  to  Dykstra’s  open-ended  deriva¬ 
tive  Hartree-Fock  (DHF)  method.^^  Others  have  now 
written  similar  TDHF  programs  which  employ  the  2n  4- 1 
rule  for  0^^  and  for  0  and  while  our  original  procedure 
permits  the  higher  polarizabilities  5,e,  •  •  • ,  to  be  readily  ob¬ 
tained  as  well.  Because  of  the  absence  of  electron  correla¬ 
tion,  the  TDHF  theory  is  not  usually  sufficient,  alone,  to 
predict  the  absolute  value  of  such  properties.  Others  have 
considered  the  more  general  formulation  of  frequency  de¬ 
pendent  hyperpolarizabilities. However,  we  observe 
that  correlated  static  values  augmented  by  the  percentage 
correction  for  the  frequency  dependence  given  by  TDHF'^ 
seems  generally  justified.  In  cases  where  the  correlated  fre¬ 
quency  dependence  is  known  from  MBPT(2),'*  or  from 
multiconfiguration  linear  response  (MCLR)  calcula¬ 
tions,^'  the  TDHF  percentage  dispersion  is  quite  close.  Re¬ 
cently,  we  reformulated  our  original  atomic-orbital-basis 
higher-order  TDHF  computational  algorithm  in  the  mo¬ 
lecular  orbital  basis  to  exploit  simplifications  in  the  Fock 
matrices  that  permit  a  more  efficient  application  of  TDHF 
theory  for  many  molecular  systems.”  The  program  also 
exploits  the  2n  rule  for  calculating  the  hyperpolarizabili¬ 
ties. 

Vibrational  corrections  are  also  sometimes  important 
in  obtaining  predictive  accuracy  for  molecular  hyperpolar¬ 
izabilities.  There  are  two  types:  the  first  simply  involves 
averaging  over  'he  zero-point  vibration;  while  the  second, 
0yi^  and  yv,h.  involves  contributions  to  the  hyperpolariz¬ 
ability  from  vibronic  intermediate  states.^ '  Both  have  been 
previously  considered  in  calculations^*'”  and  experi¬ 
mentally.” 


In  the  following,  after  a  discussion  of  the  computa¬ 
tional  methods,  we  present  results  for  the  molecules  H-, 
N.,  CO,  CO,,  C,H^.  FH,  HT).  NH,.  and  H;S  with  aug¬ 
mented,  polarizability  consistent  basis  sets  deriscd  from 
Sadlej.'*" 

COMPUTATIONAL  PROCEDURES 

The  following  results  obtained  for  static  and  frequency 
dependent  hyperpolarizabilities  are  obtained  with  the  .aci;s 
II  program  system.*'  See  Ref.  16  for  the  relevant  equations. 
SCF  static  results  are  augmented  by  correlated  r<-sults  at 
the  MBPT(2),  CCSD.**  and  CCSDfT)*'”"'  levels.  At  the 
SCF  level  all  hyperpolarizabilities  are  obtained  analytically 
from  our  TDHF  program  for  oj-0.  At  the  correlated  lev¬ 
els,  depending  upon  symmetry  for  the  field  dependent  cal¬ 
culations.  we  either  evaluate  the  dipole  moment  in  the 
presence  of  an  electric  field,  analytically,  from  which  sub¬ 
sequent  numerical  differentiation  provides  static  a,  0.  and 
y;  or  ','ke  high-order  energy  derivatives.  In  Acts  11,'  we 
also  can  compute  analytical  second  derivatives  (i.e.,  a)  at 
the  MBPTi,2)  level,'"'”  which  permits  0  and  y  to  be  ob¬ 
tained  by  subsequent  numerical  differentiation.  Field 
strengths  are  chosen  to  be  0.005  a.u.,  which  we  have  found 
acceptable  in  previous  studies.*  The  direct  computation  of 
derivatives  is  numerically  preferable  to  a  sum  over  stales 
(SOS)  method  since  no  truncation  of  states  need  be  made. 
The  interpretative  value  of  the  two-  and  three-state  mod¬ 
els'*^  that  arises  from  the  SOS  procedure  is  lost,  however. 

We  introduce  frequency  dependence  for  the  various 
optical  processes  via  TDHF  calculations  at  the  appropriate 
o).  We  analytically  compute  all  the  different  possible  pro¬ 
cesses  and  components  for  0  and  y,  without  using  Klein- 
man  symmetry, The  correlated  results  are  scaled  by  the 
TDHF  percentage  dispersion.  For  cases  where  the  corre¬ 
lated  dispersion  is  known,'*"*'  we  observe  that  for  NH^  the 
SHG-TDHF  dispersion  at  0.06  a.u.  (clos  ■  to  the  ruby  laser 
frequency  of  0.0656  a.u.)  gives  359r,  which  is  exactly  that 
for  MBPT(2)  and  17%  compared  to  23%  for  H,CO.  Sim¬ 
ilarly,  for  EOKE,  TDHF  gives  10%  compared  to  9.8%  in 
.second  order  for  NH,,  and  5.1%  compared  to  6.4%  for 
H2CO."*  For  FH,  the  dispersion  of  the  MCLR  SHG  results 
at  w  =  0.06  a.u,  also  appears  to  be  effectively  parallel  to 
those  from  TDHF."  For  y  le.ss  correlated  dispersion  infor¬ 
mation  is  available.  For  FH,  MCLR  calculations  report 
19%  for  dc-SHG  compared  to  15%  obtained  by  TDHF  ' 
The  other  errors  in  the  calculations  reported  are  such  that 
a  few  percent  difference  between  TDHF  and  the  correct 
dispersion  will  not  change  any  conclusions  regarding  our 
results.  Of  course.  TDHF  and  correlated  frequency  depen¬ 
dent  methods  will  fail  near  a  resonance  so  such  situations 
are  excluded.  All  results  reported  here  are  for  the  ruby 
laser  frequency. 

Except  where  noted,  the  molecular  geometry  is  as¬ 
sumed  to  be  at  the  experimental  equilibrium,  shown  in 
Table  I.  In  each  case  our  initial  basis  is  chosen  to  be  the 
polarizability  consistent  basts  (POL)  set  recommended  by 
Sadlej.'*’  To  report  results  for  a  large  number  of  cx.imples, 
we  choose  this  basis  set  because  it  has  motiesi  size  yet 
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TABLE  I.  Molecular  geometries  and  location  of  lone  pair  basis  func¬ 
tions ' 

CO,  1.612  26  A--  2.1944  a.u. 

H;6  C.,;  /?,, ,,  =  0.957  A  /  0|,„„  ■=  104.5-. 

Two  lone  pair  basis  sets  are  located  in  the  plane  perpen¬ 
dicular  to  the  molecular  plane  with  a  distance  of  0.883  A 
from  O  and  the  angle  between  the  lone  pair  being  1 18“ 

HE  C,.,  i«HF  =  0.917  A=  1.7328  a.u.;  /tVif  =  0,932  =  1 .7621 
a.u. 

Three  lone  pair  basis  sets  are  located  at  the  position  of  the 
bond  distance  R^f  with  terahedral  angles. 

NHj  C,,t  =  1 .012  a.  Z  ©hnh  =  106,7“. 

One  lone  pair  basis  set  is  located  in  the  molecular  axis 
with  the  bond  distance  Rmt- 
CO  C,,;  Rco=1128^2  A  =  2.I32  22!  a.u, 

N,  R^.^,..=  1.097  68  A  =  2.074  32  a.u. 

O,*;R„h=0T4O847  97  A=l.4a.u. 

H;S  C-  t  /?HS=  1-3280  A,  Z©hsh  =  ^2.2";  two  lone  pair  basis 
sets  are  t',— -•»-<  m  the  plane  perpendicular  to  the  molec¬ 
ular  plane  at  a  distance  of  1.328  A  from  5  and  the  angle 
between  the  lone  pairs  being  92.2°. 

C2H4  Rcc“1-34  A,  /  ©(^cc  — 2  ©hch”  120°;  Rhc'~108 

A. 


‘Spectroscopic  Data  Relative  to  Diatomic  Molecules  (Pergamon,  New 
York,  1970);  Tables  of  Interatomic  Distances  and  Configuration  in  Mol¬ 
ecules  and  Ions,  Special  Publication  No.  18  (The  Chemical  Society,  Lon¬ 
don,  1965). 

almost  the  same  quality  as  larger  ANO  basis  sets'*^  in  its 
ability  to  reproduce  the  molecular  dipole  polarizabilities  in 
correlated  calculations,^^  We  identify  the  POL  basis  as  ba¬ 
sis  set  1.  Basis  set  2  is  created  from  basis  set  1  by  augment¬ 
ing  the  POL  basis  with  a  set  of  d-type  Cartesian  Gaussian 
functions  (^^=0,1)  for  hydrogen.  The  exponent  is  chosen 
to  reproduce  the  hydrogen  molecule  hyperpolarizability. 
Since  the  basis  is  intended  to  be  used  in  larger  molecular 
systems,  only  one  set  of  extra  functions  is  added.  As  will  be 
seen,  these  extra  polarization  functions  are  very  important 
in  describing  the  high  polarization  of  the  hydrogen  atom. 

For  the  polar  molecules  FH,  H2O,  NH3,  and  HjS,  we 
use  an  unconventional  type  of  basis  as  basis  set  3.  This 
basis  adds  to  basis  set  2  lone-pair  functions  which  consist 
of  an  s  and  a  set  of  p  functions  (fjp=0. 1 )  located  in  the 
lone-pair  direction  (see  Table  I)  with  the  same  distance 
from  the  heavy  atom  X,  as  X-H  where  X  =  N,  O,  F,  or  S. 
Our  objectives  in  introducing  lone-pair  functions  include 
the  following; 

( 1 )  To  introduce  a  dramatic  change  compared  to  con¬ 
ventional  atom  centered  basis  sets  to  partly  assess  the  de¬ 
gree  of  basis  set  convergence; 

(2)  to  alleviate  the  necessity  of  adding  many  high  an¬ 
gular  momentum  polarization  functions  to  the  heavy  atom; 

{ 3 )  to  provide  a  description  more  .specific  to  lone  pairs, 
which  are  generally  less  well  described  in  more  conven¬ 
tional  quantum  chemical  basis  sets.  Comparisons  with  ex¬ 
tensive  atom  centered  basis  sets  are  considered  in  the  Dis¬ 
cussion. 

RESULTS  AND  DISCUSSION 
COj 

For  the  nonpolar  CO2  molecule,  results  are  shown  in 
Table  11.  The  first  observation  is  that  the  SCF  components 


of  the  polarizability  in  the  PGL  basis  are  quilt*  similar  lu 
the  large  basis  (5i'4p3d)  f  )  results  of  Amos  1.7  al.^'  who 
find  a„=23.7  and  o,,.-il.9.  Other  extensive  basis 
(te4.p4d!/  )  results  of  Maroulis  anti  Thakkar  give  23.7 
and  12.0.'*'*’*  Visser  and  Wormcr  in  a  (  6v5p3d2/)  basis 
obtain  23.6  and  I  l.d,'*'’  as  well.  Correlation  has  a  significant 
effect  of  about  10%  giving  a  —  17.6  fivr  CCSD(T),  which 
agrees  well  with  experiment  and  with  the  SDQ-MBPT(4) 
value.'***’’  SDQ(4)  is  the  fourth-order  approximation  to 
CeSD,  which  offers  the  more  appropriate  comparison, 
however.  As  usual,  the  anisotropy  is  less  accurate,  but  our 
CCSD(T)  value  of  14.5  (or  CC'tf)  ^vf  14.4)  compares  quite 
well  with  the  SDQ(4)  result  of  U.J."***” 

The  static  SCF  second  hyperpolarizability  y  is  800 
wh;,h  is  about  60%  of  the  dc-SHG  experimental 
value  of  1332  a.u..  MBPT(2)  improves  the  result  at  the 
static  limit  to  1170  a.u.,  i.e.,  a  ~30%  improvement.  The 
effect  of  infinite-order  correlation  as  measured  by  CCSD 
actually  reduces  the  correlation  effect  to  give  1070  a.u. 
Introducing  TDHF  dispersion  increases  y  by  31%  (or 
accounts  for  13%  of  the  experimental  value)  for  THG  and 
13%  (or  8%  of  experiment)  for  dc-SHG.  Scaling  by  the 


TABLE  U.  Components  of  the  polarizability  and  second  hyperpolanz- 
abilities  of  the  COj  molecule  and  results  for  various  processes  (in  a  u,).* 


SCF/ 

TDHF" 

MBPT{2)* 

C..SD* 

CCSD(T)» 

EXP 

11.8101 

12.8638 

12  6327 

12,7908 

“z. 

7.3.8862 

27  8482 

27,0631 

27,2984 

a' 

15.835 

17,859 

17.443 

17,627 

17  50' 

Aa’’ 

12076 

14.984 

14430 

14.508 

13,83' 

700  0 

1010  0 

9300 

900.0 

810.0 

1120,0 

9000 

10.100 

Yma 

330.0 

500,0 

4900 

520,0 

n  '■ 

800.0 

11700 

10700 

1150,0 

THG" 

1046.0 

1500.0 

1400.0 

1500.0 

1860i:270‘ 

dc-SHG' 

910  0 

1300,0 

1200.0 

1300.0 

1332  7  IS"" 

IDRl* 

871  0 

1.100.0 

1200.0 

1200  0 

EOKE" 

8330 

1200.0 

11000 

1200.0 

“The  2  axis  is  the  principal  axis  of  the  molecule. 

'’In  this  and  subsequent  tables,  frequency  dependent  polarizabilities  are 
evaluated  at  w=0  0656  a  u.  =  6443  A. 

‘iT  =  -j(a„.f  a.,  +  a„). 

=  ((a„-a,,)%  («„ -  n„)%  7 

~  iYu/;+  Einstein  summation  convention. 

A-  Spackman,  J.  Phys.  Chem.  93,  7594  (1989).  obtains  these  values 
from  analysis  of  several  dilferenl  experiments. 

*Corrclatcd  frequency  dependent  values  are  estimated  by  y 
r-{y,j  (ORR^^I  ‘n  this  and  subsequent  tables 

’’THG  third  harmonic  generation  y(  - 
'dc-SHG:  dc-indiiced  second  harmonic  generation,  y( 

'IDRI;  intensity  dependent  refractive  index.  y(  ^  <>}) 

’“EOKE:  electro-optic  Kerr  effect.  y(  dc-OR  de- induced  op¬ 
tical  rectification,  y(0;fu.  is  numencaily  the  same  as  HOKE  after 

averaging 

Reference  5.  Since  y  /6  conversion  factor  is  ^  v  5  0.i  /  \  10  *’  esu/ 
mclcculc/a  u.  7  8  .^95  10  esu/molcculc/a  u 

"'Reference  4 


J.  Chem  Phys.,  Vol.  90,  No.  4,  15  February  1993 


K  Sekino  and  R  Bartlett:  Molecular  hyperpoianzabihiies 


3025 


TABLE  III.  Components  of  the  polarizability  and  second  hyperjKilari/- 
abilities  of  the  N-  molecule  and  results  for  various  processes  (m  a  u.  j  ' 


SCF/ 

TDHF 

MBPT(2) 

CCSD 

CCSDtT) 

EXP 

9,6310 

>),9750 

10.0002 

10  1290 

15.0033 

14.4138 

14.8266 

14  9975 

a 

11.422 

11.455 

11.609 

11.752 

11.8" 

Aa 

5.372 

4.439 

4.826 

4.869 

4  67' 

y  xjtxj 

660.0 

7800 

8100 

860.0 

790.0 

1220.0 

1220.0 

1290.0 

VxAi; 

270.0 

340  0 

350.0 

370.0 

T; 

730.0 

930.0 

950.0 

1010.0 

THG' 

937.0 

1200.0 

1200.0 

1300.0 

1 270.(7’ 

dc-SHG 

822.0 

1100.0 

1100.0 

1100,0 

1030 ± 12“ 

IDRI 

788,0 

1000.0 

1000.0 

1100,0 

EOKE 

756.0 

1000.0 

1000.0 

1 100.0 

“The  z  axis  is  the  principal  axis  of  the  molecule.  )-  .5,  Aa  are  defined  in 
the  footnotes  of  Table  II 

‘’A,  D.  Buckingham,  M.  P.  Bogaard,  D.  A.  Dunmur,  C.  P  Hobbs,  and  B- 
J.  Orr,  Trans.  Faraday  Soc.  66.  1548  (1970). 

'All  processes  are  defined  in  the  footnote  of  Table  11. 

‘’Reference  4. 


TDHF  percentage  correction  to  the  correlated  static  re¬ 
sults  {MBPT(2)  and  CCSD]  for  dc-SHG  gives  I30()  and 
1200  a.u.,  respectively,  while  triples  as  measured  by 
CCSD(T)  give  1300  a.u.,  and  are  in  good  agreement  with 
the  experimental  value  of  1332  ±15  a.u.  The  theoretical 
estimates  for  THG  of  15(X),  1400,  and  15(X)  a.u.  are  lower 
than  the  experimental  mean  value,  I860 ±270  a.u.  but 
nearly  fall  within  the  rather  large  error  bars.  The  other 
correlated  [SDQ(4)]  study  of  the  hyperpolarizability  of 
gives  for  static  1 197  a.u.  compared  to  our  1 150 
a.u.,  showing  excellent  agreement.  Our  results  suggest  the 
triples  would  likely  increase  that  value  by  7%.  For  CO2, 
the  Tvib  ^Iso  been  determined  to  be  —86.83  at  <d  =  0.07 
a.u.  for  dc-SHG  and  +2.0  a.u.  for  THG.^'' 

Nj 

The  Nj  results  are  shown  in  Table  III.  The  behavior  is 
quite  similar  to  that  for  CO2.  At  the  SCF  [i.e.,  TDHFfiu 
=  0)]  level  we  obtain  5=11.42  compared  to  11.40  ob¬ 
tained  in  a  (6s5p3cl2  f)  basis'**'^'’’  and  11.57^*^’  in  a 
{6s^p3d\ f  )  basis.  Our  SCF  Aa  of  5.37  compares  reason¬ 
ably  to  5. 1 1  and  5. 17  from  those  calculations.  As  we  intro¬ 
duce  correlation  our  a  increases  to  be  in  good  agreement 
with  experiment,  while  Aa  appropriately  decreases.  Other 
SDT-MBPT  results  behave  differently  increasing  Aa  to 
5.48.‘**  It  should  be  noted,  however,  that  this  particular 
MBPT  approach  presents  a  diagrammatic  evaluation  of 
second-order  property  diagrams  (see  Ref.  47  for  a  discus¬ 
sion)  and  is  not  equivalent  to  the  results  from  a  finite-field 
energy  derivative  approach  which  would  introduce  a  larger 
selection  of  correlation  contributions.  In  particular,  using 
the  latter,  SDQ(4)  gives  4.65.'**“^ 

The  static  yy  SCF  value  is  730  a.u.  Correlation  effects 
increase  this  by  up  to  40%,  with  triples  being  6%  of  the 
value,  compared  to  7%  for  COy.  In  the  CO2  case,  CCSD 
reduces  the  MBPT(2)  y.j  instead  of  enhancing  it.  The  950 


a.u,  CCSD  value  is  in  reasonable  agreenicm  vviih  the 
SDQ(4)  value  of  830.’'  "  l  ike  correlation.  1  DHF  disper¬ 
sion  efiects  show  a  similar  incremen!  t'f  3%  to  3!)'T  of  l.he 
static  values,  m  ihe  uvual  order  IMG  •  dc-SHe'  -  IDRI 
>  EOKE.  Introducing  correlation  and  dispersion,  reason¬ 
able  agreement  with  THG  experiment  is  obtained,  al¬ 
though  unlike  CO,  there  is  somewhat  more  error  for  dc- 
SHG.  Results  for  optical  properties  at  several  different 
frequencies  are  presented  elsewhere.'" 

H2 

The  polarizabilities  and  the  static  second  hyperpolar- 
izabilities  of  are  presented  in  Table  IV.  As  our  results 
are  only  for  the  electronic  part,  we  compare  with  the  exact 
electronic  results  taken  from  the  very  accurate  explicitly 
correlated  calculation  of  Bishop  et  where  the  addi¬ 
tional  vibrational  effects  are  also  presented.  Our  relatively 
poor  agreement  with  the  exact  electronic  results  requires 
some  basis  set  extension.  For  example,  the  calculated  ratios 
Yizz/Yxxxx  set  1  are  more  than  3.  By  adding  the 

diffuse  polarization  (d-function  exponent  0.1).  de- 
crca,ses,  while  increases  and  the  ratio  of  the  two  com¬ 
ponents  becomes  dose  to  one.  The  y,,..  component  seems 
to  be  less  sensitive  to  the  choice  of  basis  set.  The  same 
tendency  in  the  improvement  of  the  hyperpolarizability 
components  was  observed  in  a  previous  study.'"  Similarly, 
even  the  polarizability  and  anisotropy  show  improvement 
from  basis  set  I  to  basis  set  2.  Basis  set  2.  despite  its  rela¬ 
tively  modest  size,  provides  the  parallel  component  y 
=  610  that  is  comparable  to  the  exact  theoretical  calcula¬ 
tion  at  the  static  limit,  613.2.^**  For  dc-SHG  at  /I  =694  3 
nm  TDHF  dispersion  correction  increases  our  static  value 
to  706.4  a.u.  The  corresponding  value  calculated  from  the 
exact  components  is  721.2  a.u.  Thus  the  percentage  cor¬ 
rection  for  the  dispersion  effect  of  15.8%  by  TDHF  is 
slightly  lower  than  that  obtained  by  the  exact  calculation 
(17.6%).  This  seems  to  be  typical  of  TDHF  as  a  measure 
of  dispersion,  as  it  usually  obtains  the  correct  effect  at 
standard  frequencies  to  within  a  couple  of  percent,  barring 
pathological  cases.  TDHF  or  even  any  correlated  treat¬ 
ment  of  dispersion,  will  certainly  fail  near  a  resonance. 
TDHF=RPA  is  also  known  to  generally  provide  reason¬ 
able  singlet  excitation  energies,  which  is  necessary  if  it  is  to 
give  a  good  dispersion  estimate. 

Table  IV  shows  the  correlated  for  the  various  pro¬ 
cesses  after  correction  for  the  dispersion.  The  correlated 
results  such  as  MBPT(2)  and  CCSD  (full  Cl  for  H,)  are 
improved  and  close  to  the  exact  results.  After  averaging 
over  the  vibrational  zero-point  wave  function  the  exact  v 
becomes  682.7  a.u.  at  the  static  limit  and  811.6  a.u.  at 
A  =694,3  nm  for  dc-SHG.  This  is  a  part  of  the  vibrational 
effect.s.  The  other  so-called  “pure”  vibrational  contribution 

(i.e.,  the  part  that  involves  scattering  through  vibronic 
states')  to  Pji  is  a  large  t  183.7  a.u.  at  the  static  limit.'  ''  " 
while  the  more  appropriate  frequency  dependent  value  is  a 
smaller  -  13.5  a.u.  at  A  694.3  nm.  This  results  in  a  y  of 
866.4  and  798.1  a.u..  respeclively.'"  As  can  be  easily  seen 
from  the  actual  calculations'"'  no  percentage  dispersion 
correction  can  be  applied  for  The  importance  of  ;y,,,  is 
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TABLE  !V.  CumfKineni<i  the  polanzability  and  second  hypcrp<>lariz- 
abiltties'  of  the  H<  molecule  and  results  for  various  processes  (in  a  u  ) 


SCF/ 

TDHF 

.MBPT(2) 

CCSD 

Exact” 

EXP 

4.1940 

4  .3655 

4,3423 

4  hl99 

•»608! 

4,.'«fi9 

4,5786 

6.5481 

6.5693 

6.5640 

6  4453 

6.4846 

6.4929 

6,3875 

a 

5,!12 

5.100 

<082 

5.228 

5.234 

5.223 

5,1816 

5.43' 

Aa 

2.154 

2.204 

2.222 

t.825 

1.877 

1.906 

2.0371 

2.035'' 

240  0 

2600 

270.0 

530.0 

580,0 

590.0 

575.9 

760.0 

800.0 

810.0 

560.0 

620.0 

630.0 

682.5 

Y-XXZZ 

220.0 

230.0 

230,0 

200.0 

210.0 

210.0 

211,9 

y-i 

450.0 

480.0 

490.0 

550.0 

600.0 

610.0 

613  2 

THG 

631.0 

670.tf 

680,0' 

747.0 

810.0 

830,0 

858.0' 

950  ±  140' 

dc-SHG 

SlQ.O 

560.0 

570,0 

639.0 

690.0 

710.0 

721.0 

801  ±4.5' 

IDRI 

505.0 

530.0 

540.0 

607.0 

66C.0 

670,0 

EOKE 

478-0 

510-0 

520.0 

578.0 

630.0 

640.0 

6460 

560  ±60' 

*The  upper  and  lower  numbers  are  obtained  by  the  calculations  using 
basis  set  1;  (3s2p]  and  basis  set  2;  [3s2pld]  (d  function  on  hydrogen  has 
exponent  0.1). 

*The  numbers  are  evaluated  by  the  formulas  in  the  footnotes  c-e  in  Table 
n  using  the  pure  electronic  values  from  Ref.  48. 

'A.  C.  Newell  and  R.  C,  Baird,  J.  Appl.  Phys.  3«,  3751  (  1965). 

'’K.  B.  MacAdam  and  N.  F.  Ramsay,  Phys.  Rev.  A  6,  898  (1972). 

'The  experimental  values  are  from  D.  P.  Shelton,  Phys.  Rev.  A  42,  2578 
( 1990),  and  include  vibrational  effects  in  addition  to  the  electronic  con¬ 
tributions.  See  the  text. 

^See  footnote  g.  Table  II. 

process  dependent,  but  it  is  quite  small  for  the  frequencies 
and  processes  we  study  here.  The  rather  large  vibrational 
correction  of  — 13.5  a.u.  in  this  case  is  peculiar  to  H2.  Of 
course,  vibrational  effects  in  all  H-containing  molecules 
would  be  expected  to  be  comparatively  significant.  Allow¬ 
ing  for  corrections,  the  computed  theoretical  numbers  ob¬ 
tained  in  basis  2  are  reasonably  close  to  the  exact  values; 
thus  we  decided  to  use  this  H-basis  set  for  a  systematic 
study  of  hyperpolarizabilities  of  larger  H-containing  mol¬ 
ecules. 

C2H4 

Ethylene  is  the  first  of  the  linear  polyene  series  whose 
nonlinear  optical  behavior  provides  a  reference  point  for 
the  development  of  optical  polymer  materials.  The  compo¬ 
nents  of  the  polarizability,  a,  and  hyperjx)larizability,  y, 
for  C2H4  are  shown  in  Table  V.  The  results  employ  basis  2, 
that  is  the  POL  basis  for  C  and  H,  with  the  latter  aug¬ 
mented  by  the  set  of  diffuse  d  functions  found  to  be  impor¬ 
tant  in  H2.  The  extra  Cartesian  d  functions  introduce  a 
linear  dependence  that  requires  removal  to  get  adequately 
converged  results.  The  agreement  with  experiment  for  a 
and  Aa  is  satisfactory,  particularly  allowing  for  slight  dif- 


TABLE  V.  Compoaenls  ot  the  {■HtlarizabilitN  and  sccoud  fivpcff>«dunz 
ubililies  of  CjH4  uHti  results  for  varuius  prt>cesscs  sin  a  «  ]  ^ 


SCF/ 

rDHi 

MBPTC) 

CCS1.5 

CeSDi.  1  ) 

EXP 

.16  37.1 

.14,0.14 

.13  814 

3.1  854 

24622 

25  141 

24.9.14 

25  148 

2.V064 

22  85.1 

22  X)}} 

22  2-.W 

a 

28  02 

27  41 

27  08 

2h  7(f' 

Aa 

12.60 

10  17 

1065 

10  47 

10  98' 

r,... 

3  3(X)  0 

6  1(X3.0 

5  4<8V0 

5  tXX)  0 

/ 

2  800.0 

3  5(X),0 

.1  2(X)  0 

4(X)  0 

Yidii 

1 1  900  0 

1 1  800  0 

9  800  0 

10  2000 

Yix:t 

3  100.0 

1  5(X).0 

2  900  0 

1  (XX)  0 

2  500.0 

2  600  0 

2  2(X>0 

2  ItX)  0 

r,.., 

1  6000 

2  (XX).0 

1  800  0 

1  9(X)  0 

r 

6  500  0 

7  500.0 

6  4(X)  0 

6  7(X)  0 

THG 

17  500  0 

20  300  0 

17  3(X),0 

18  200  0 

dc-SHG 

9  900  0 

1 1  400  0 

9  7000 

10  2tX)  0 

9030  t  2tXE 

IDRI 

8  500  0 

9  8fX).0 

8  400  0 

8  8(X).0 

EOKE 

7  400.0 

8  600.0 

7  300.0 

7  7(X)  0 

*The  X  axis  i.s  the  pnnctpal  axis  of  the  molecule,  the  y  axis  is  in  the  plane 
and  the  z  axis  is  out  of  the  plane  of  the  molecule  Basis  set  2 
•T.  K.  Bose  and  R.  H  Cole.  J.  Chem  Phys.  54.  3829  ( 1971  ) 

'Reference  5. 


ferences  from  our  vibrationless  results  and  the  experimen¬ 
tal  values  that  are  obtained  including  zero-point  effects. 

Introducing  the  frequency  dependence  (Table  V),  the 
TDHF  result  is  in  reasonable  agreement  with  the  dc-SHG 
experimental  value.  However,  MBPT(2)  increases  the 
value  significantly.  As  occurs  in  CO2,  the  infinite-order 
CeSD  method  decreases  the  second-order  result,  and  in 
this  case,  even  to  a  value  less  than  the  TDHF  result.  Triple 
excitations  increase  the  value  slightly  leaving  the  final  dc- 
SHG  result  in  excess  of  experiment  by  13.4%.  The  corre¬ 
lation  effects  in  ethylene  are  the  smallest  we  have  seen. 
They  are  larger  in  fzafls-butadiene,'*  '’  the  next  member  of 
the  polyene  series.  The  other  optical  processes  show  similar 
behavior.  This  is  the  one  example  we  have  seen  where 
MBPT(2)  significantly  overestimates  the  correlation  cor¬ 
rections  to  the  hyperpolarizability.  Larger  basis  studies  are 
warranted. 

CO 

The  CO  molecule  has  a  jS  hyperpolarizability  along 
with  y.  These  are  shown  in  Tables  VI  and  VII.  Along  with 
/?,  we  also  report  which,  because  of  its  small  value 
changes  sign  from  the  SCF  value.  For  CCSD  provides 
excellent  agreement  with  the  experimental  value  although 
including  triple  excitations  via  CCSD(T)  hurts  the  agree¬ 
ment.  Slight  basis  set  changes  will  affect  the  dipole  moment 
in  the  third  decimal.  For  the  polarizability,  many  extensive 
basis  set  calculations  have  been  reported,  and  are  well  sum¬ 
marized  by  Rijks  and  Wormer.'*'*  Our  result  for  a  is  in 
excellent  agreement  with  the  vibrationally  corrected  exper¬ 
imental  result  quoted  by  Werner  and  Meyer.'"  Lhc  anisot¬ 
ropy,  as  usual,  is  somewhat  poorer  but  consistent  with 
Wernei  and  Meyer’s  CEPA  value  of  3.91.  Sunil  and  Jor¬ 
dan”  get  3.79  using  CCD. 
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TABLE  VI.  Dipole  moment  and  components  of  the  first  hyperpolarie- 
abilities  of  the  CO  molecule  and  results  for  various  processes  (in  a  u  ) 


SCF.TDHF 

MBPT(2) 

CCSD 

CCSD(T) 

EXP 

p.- 

-0.0989 

-f  0.1 208 

+0.0421 

f0.n581 

+  0  044'- 

{-)28,2 

25.6 

26.1 

26.0 

ft.. 

(-)3.5 

6.0 

6.1 

6.6 

0 

(  -)21.1 

22.6 

23.0 

23.5 

SHG-* 

24.1 

25.9 

26.4 

27,0 

29,9  +  3  2' 

OR  or 
EOPE' 

21.9 

23.6 

24.1 

24  6 

“The  z  axis  is  the  principal  axis  of  the  molecule. 

''/3ii  =(.0,ji+0j,j+Pjj,)/i  with  summation  convention. 

‘^Reference  50,  and  references  therein. 

■^SHG:  second  hannonic  generation,  (i( 

'Reference  4.  Since  Y'’f’  =  |3s  /2  conversion  factor  is  |x8.639x  10  ” 
=4.3195  X  10“’*  esu/molecule/a.u. 

'OR;  optical  rectification.  i9t0:&i,  —  in).  HOPE:  electro-optic  Pockets  effect. 
0(  —  t!>;6j,0).  After  averaging  these  two  are  numerically  equal. 

Discussing  the  yj|  hyperpolarizability  first  we  see  some 
similarities  with  CO2  and  Nj.  The  static  Y\i  value  is  in¬ 
creased  via  electron  correlation  by  up  to  57%.  Like  COj, 
and  unlike  Nj,  the  infinite-order  correlation  effects  mea¬ 
sured  by  CeSD  slightly  reduce  the  correlation  correction, 
while  triples  further  enhance  it.  The  7%  effect  of  triples  is 
consistent  with  that  observed  for  N2  and  CO2.  TDHF  dis¬ 
persion  varies  from  6%  for  EOKE  to  47%  for  THG.  Com¬ 
bining  dispersion  and  correlation,  the  dc-SHG  result  is 
slightly  too  large  compared  to  experiment  as  is  that  for  N2. 

The  first  (/?)  hyperpolarizability  with  dispersion  al¬ 
lows  the  evaluation  of  second  harmonic  generation.  Pi 
—  2io-,(o.ci})  (SHG);  optical  rectification,  P(0;a,—<o) 
(OR);  or  the  electro-optic  Pockels  effect,  Pi—(o;o>,0) 
(EOPE).  For  )3||  —iPm-YPjij+Pjjp/^  (summation  con¬ 
vention)  the  latter  two  are  numerically  equivalent.  Unlike 
y,  P  has  a  sign,  and  is  defined  to  be  positive  when  the 

TABLE  VIL  Components  of  the  polarizability  and  second  hyperpolariz¬ 
abilities  of  the  CO  molecule  and  results  for  various  processes  (in  a.u.).“ 

SCF/ 


TDHF 

MBPT(2) 

CCSD 

CCSD(T) 

EXP 

11.1068 

11.7888 

11.5997 

1 1,7332 

«« 

14.4752 

15.7020 

15.5647 

15,6522 

a 

12.23 

13.09 

12.92 

13.04 

U.og" 

Aa 

3.37 

3.91 

3.97 

3.92 

3.65' 

y xxxx 

920.0 

1380.0 

1360.0 

1470.0 

7zZ22 

1200.0 

1740.0 

1740.0 

1880.0 

Yxxzx 

360.0 

520.0 

510,0 

540.0 

rn 

1020.0 

1500.0 

1480.0 

1590,0 

THG'* 

1484.0 

2200.0 

2200.0 

2300.0 

dc-SHO 

1211.0 

1800,0 

1800.0 

1900.0 

1720  +  48' 

IDRI 

1138,0 

1700.0 

1700.0 

1800.0 

EOKE  or 

1071.0 

1600.0 

1600.0 

1700.0 

OR 

The  2  axis  is  the  principal  axis  of  the  molecule. 

"Reference  50,  includes  zero-point  ground  state  vibrational  correction. 
‘G.  H  F.  Diercksen,  B  O.  Roos,  and  A.  1  Sadlej,  Chem.  Phys,  59,  29 
(19811. 

‘’All  quantities  are  defined  in  the  footnote  m  Table  II,  and  correlated 
results  incorporate  the  percentage  TDHF  dispersion  corrccuon. 
'Reference  4. 


product  P  • /i  is  positive,  that  is  when  it  has  the  same 
direction  as  the  dipole  moment  The  evaluation  of  P  is 
simpler  than  y  in  the  sense  that  it  is  a  lower-order  (third- 
order)  quantity  in  the  electric  field,  and  only  requires  a 
knowledge  of  ilie  first-order  perturbed  wave  funciitin  if  the 
usual  2n+  I  rule  of  perturbation  theory  is  used.  However, 
it  is  more  complicated  in  the  sense  that  it  depends  upon 
delicate,  relative  differences  in  the  molecule’s  charge  den¬ 
sity,  causing  it  to  be  more  sensitive  to  different  basis  set 
characteristics  than  is  y. 

Since  the  dipole  moment  of  CO  is  known  to  have  the 
wrong  sign  at  the  SCF  level,  which  is  corrected  once  an 
adequate  amount  of  electron  correlation  is  inciuded,  the 
SCF  value  of  p,  must  have  the  opposite  sign  to  /i^ci ,  but 
the  same  sign  once  correlation  is  introduced.  This  is  what 
we  see  as  indicated  by  the  (  —  )  sign  in  the  SCF  column  of 
Table  VI.  In  other  words  the  sign  of  is  correct  even 
though  that  of  not. 

The  correlation  corrections  to  P  show  a  more  mtxlesi 
increase  than  was  observed  for  ,  being  1 1%  of  the  static 
value  at  the  CCSD(T)  level.  Dispersion  also  increases  the 
result  by  4.4%  for  OR  and  up  to  14%  for  SHG.  Combin¬ 
ing  the  dispersion  and  correlation  effects,  for  this  relatively 
nonpolar  molecule,  the  agreement  with  experiment  for 
SHG  is  quite  good.  The  basis  set  is  certainly  not  converged, 
however,  and  extensions  might  change  the  degree  of  agree¬ 
ment. 

HF 

In  Tables  VIII  and  IX,  we  present  the  dipole  moment, 
polarizabilities,  and  hyperpolarizabilities  of  the  FH  mole¬ 
cule.  Values  are  reported  for  each  component  in  three  dif¬ 
ferent  basis  sets.  The  first  is  for  Sadlej’s  standard  POL 
basis^*  (basis  set  1),  which  appears  to  have  adequately 
described  the  yj  hyperpolarizability  for  N2,  CO,  and  COi. 
The  second  adds  the  diffuse  d  polarization  function  on 
hydrogen  (basis  set  2),  which  was  important  m  H2.  To 
potentially  introduce  a  dramatic  change  compared  to  the 
nine  rather  extensive  is,p,d,f  )  conventional  basis  sets  con¬ 
sidered  by  us  previously^'*  and  other  even  more  extended, 
atom  centered  basis  sets  considered  by  others, we 
choose  to  introduce  three  sets  of  lone-pair  functions  ( one  y 
and  a  set  of p  functions  for  each  lone  pair)  to  basis  set  2  to 
define  (basis  set  3).  As  may  be  observed  in  Table  VIII.  for 
the  dipole  moment  the  improvement  is  negligible  for  both 
extensions,  but  for  /3|j  both  effects  are  significant.  At  the 
SCF  level  each  extension  has  comparable  importance,  but 
the  synergistic  effect  with  correlation  is  even  more  dra¬ 
matic.  The  first  extension  causes  about  an  18%  CCSD 
increase  while  the  second  adds  another  8%.  The  net  effect 
is  a  24%  increase  in  \Pf.  [  from  that  obtained  with  just  the 
POL  basis.  Once  again,  the  effect  of  triple  excitations  as 
measured  by  CCSD(T)  .slightly  hurts  the  (probably  fortu¬ 
itous)  agreement  with  experiment  of  CCSD  for  the  dipole 
moment.  We  should  always  recognize,  however,  that  the 
perturbative  (T)  correction,  though  generally  good,  is  tiol 
nece.ssarily  the  correct  measure  of  triples,  as  it  contains 
only  one  of  five  fifth-order  contributions  involving  triples 
and  connected  quadruples  that  have  about  the  same  mag- 
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TABLE  Vlll.  Dipole  moment  and  components  of  the  first  hyperpolarir- 
abilities  of  the  HF  molecule  and  results  for  various  priKesses  (in  a  u  )  ' 


SCF/TDHF 

MBPT(2) 

CCSD 

CCSD(T) 

EXP 

(0 

0.7566“ 

0.7079 

0.7065 

0.7570" 

0.7082 

0.7070 

O.TSM" 

0.7069 

0.7062 

0  7008 

0.708'' 

-8.14 

-8.80 

-8.91 

-8-15 

-8.96 

-9,02 

-8.40 

-9.10 

-  9.24 

-9.62 

-0.26 

-0.23 

-0.08 

-0.02 

-1.03 

-0.79 

-0.28 

-  i.25 

-1.03 

-  1,27 

e. 

-4.58 

-5.56 

-5  44 

-4.87 

-6.61 

-6.36 

-5.38 

-6.96 

-6.78 

-7.30 

SHG 

-5.024 

-6.1 

-6.0 

-5.359 

-7.3 

-7.0 

-5.853 

-7.6 

-7.4 

-8.0 

-10  9±0  95' 

OR  or 

-4.71! 

-5.7 

-5.6 

EOPE 

-5.018 

-6.8 

-6.6 

-5.522 

-7.2 

-7,0 

-7.5 

“Basis  set  1;  [5s3p2r//3j2/)]. 

"Basis  set  2;  [isipld/is2p\d],  the  d  function  on  hydrogen  has  an  expo¬ 
nent  of  0.1, 

"Basis  set  3;  Basis  set  2  plus  lone-pair  functions  of  s  and  p  types  with 
exponent=0. 1  at  0,916  946  A  from  fluorine 
‘'Reference  50. 

'Reference  56. 

nitude,'*^  and  even  higher-order  contributions  will  have  a 
role  in  definitive  answers.  It  is  primarily  an  empirical  ob¬ 
servation  that  CCSD(T)  generally  woiKs  well  for  certain 
classes  of  problems,’^  but  it  can  fail  in  others.^  ” 

From  Table  IX,  a  similarly  modest  change  with  basis 
in  a  and  Aa  is  observed.  All  three  basis  sets  give  effectively 
perfect  agre’ment  with  other  large  basis  SCF  results  for  a 
and  for  Aa,'*'*’”  and  give  rather  consistent  results  at  all 
correlated  levels.  The  anisotropy  of  1.10  does  not  compare 
well  to  experiment,  or  with  some  other  correlated  calcula¬ 
tions  that  give  values  that  include  1.16  (MCSCF^’),  1.27 
(CEPA^°),  and  1.37  from  our  older  SDQ-MBPT(4)  re¬ 
sults.*  However,  the  CCSD,  CCSD-(-T(CCSD),  and 
CCSD(T)  results  of  Chong  and  Langhoff'^  with  an  exten¬ 
sive  ANO  basis  give  1.18,  1.15,  and  1.16,  respectively,  and 
1.18  for  MCPF,  which  are  closer  to  the  current  value. 

The  change  with  basis  is  much  more  pronounced  in 
Yii  ,  particularly  for  the  component,  which  increases 
by  about  100  a.u.  in  SCF  and  about  200  a.u.  in  correlated 
calculations.  The  results  with  basis  set  3  for  )3j|  and  yij  are 
comparable  or  even  larger  in  magnitude  than  those  with 
the  largest  basis  sets  presented  in  our  previous  work,*  *  but 
there  is  still  considerable  difference  from  experiment.’*  A 
detailed  assessment  of  the  present  results  and  comparisons 
with  other  results,  including  new  atom  centered  basis  re¬ 
sults,  will  be  considered  in  the  Discussion  section. 

HjO 

Previously,  we  presented  correlated  SDQ-MBPT(4) 
results  for  H2O  in  several  different  basis  sets,^  including  a 
(_6s5p4d/4s2p)  basis  ( D).  In  Tables  X  and  XI,  we  summa¬ 
rize  the  dipole  moment,  polarizabilities,  and  hyperpolariz- 


FABLE  I.V.  Compimcnts  iif  the  ptibri/ubilily  ami  vccoiid  hyjvcrpolati/- 
ahilitiex  of  ihc  HF  niolvcult-  and  rcsulis  for  variouv  pnivcsscs  >  in  1 11  j 


SCF/TDHF 

M  BP  1(2) 

CCSD 

cesDt  r ) 

EXP 

4.4567 

5  307 1 

5  1656 

4.4764 

5,3508 

5  1996 

4.5103 

5.3975 

5  2.190 

5.,VA9}s 

5  08- 

a.. 

5-7416 

6.4072 

6.3268 

5.7758 

6.4572 

6.  .5681 

5,7577 

64326 

6.3440 

6  4378 

bAir 

a 

4  89 

5  67 

5  55 

4.91 

5  72 

5  59 

4.93 

5.74 

5,61 

5.7) 

5.52* 

Aa 

1.28 

1.10 

1,16 

1  30 

Ml 

1  17 

1.25 

1,04 

111 

1  10 

1,32' 

250.0 

430,0 

410  0 

280,0 

500.0 

470.0 

350.0 

640.0 

600  0 

650  0 

r-v.^v 

2500 

3400 

3200 

280,0 

380.0 

3700 

270.0 

380.0 

360.0 

390,0 

700 

120  0 

110,0 

900 

1600 

150,0 

100.0 

170  0 

160  0 

180.0 

Y- 

2400 

3900 

370.0 

280.0 

470.0 

440,0 

320.0 

560.0 

520.0 

5600 

THG 

303,0 

490  0 

470.0 

347,0 

590.0 

5500 

4040 

700.0 

660.0 

710,0 

dc-SHO 

270.0 

440.0 

420.0 

830  *  120" 

308.0 

520.0 

490.0 

359,0 

620.0 

580.0 

6300 

IDRI 

260.0 

420.0 

400.0 

297.0 

500.0 

4700 

345.0 

600.0 

560.0 

610.0 

EOKE  or 

251,0 

410.0 

390.0 

dcOR 

286,0 

480.0 

460.0 

333.0 

580.0 

540.0 

5900 

*G.  H.  F.  Diercksen  and  A.  J.  Sadlej,  J.  Chem,  Phys  75,  1253  (1981). 
Experimental  values  with  zero-point  vibrational  corrections. 

"Reference  54. 


abilities  of  H2O  in  the  standard  POL  basis  1  and  basis  3. 
The  two  sets  of  lone-pair  functions  to  H,0  (basis  set  3)  are 
located  at  the  normal  OH  distance  but  directed  symmetri¬ 
cally  away  to  help  in  the  description  of  the  lone  pairs  on  O 
(see  Table  1  for  lone-pair  basis  orientation).  We  also  in¬ 
vestigated  locating  the  lone-pair  functions  at  the  covalent 
radius  of  O  but  there  were  no  important  changes.  Com¬ 
pared  to  the  older  basis  D  results,  there  is  little  change  in 
SCF  or  a  with  either  basis.  Here  CCSD(T)  gives  a 
slightly  better  than  does  CCSD.  The  very  small  anisot¬ 
ropy,  Aa,  is  off  from  experiment  by  almost  a  factor  of  2. 
The  MBPT(4)  anisotropy  from  Maroulis”*"’  is  a  similar 
0.32,  while  that  of  Diercksen  and  Sadlej  is  0.45.’' 

The  components  which  have  .small  amplitudes  .such  as 
of  HjO  are  very  sensitive  to  the  choice  of  basis  sets, 
but,  typically,  /fp  itself  is  noi  as  sensitive  as  some  of  its 
components.  The  SCF  /3i|  was  -  9.2  a.u.  in  basis  D.  Sim¬ 
ilarly,  SDt^-MBPT(4)  results  for  basis  D  gave  -  13.7  for 
^i|  falling  between  the  current  POL  values  and  those  for 
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TABLE  X.  Dipole  moment  and  components  of  the  first  hyperprdarizabil- 
ities  of  the  H2O  molecule  and  results  for  various  processes'  (in  a  u  ). 


SCF/TDHF 

MBPT(2) 

CeSD 

CCSDtT) 

EXP 

p.- 

0.7804*’ 

0.7297 

0.7287 

0,7793'' 

0.7272 

0.7273 

0.7196 

0.721'* 

-4.6 

-8.7 

-7,3 

-7.1 

-13.7 

-12.0 

-13.7 

0^ 

-0.6 

-2.4 

-2.0 

-1,4 

-5.9 

-5.2 

-6.2 

&zyy 

-9.7 

-10.6 

-lO.g 

-9.4 

-9.4 

-9.9 

-10  2 

0-:l 

-8.3 

-13.0 

-12.1 

-10.8 

-17.5 

-16.2 

-18,0 

SHG 

-9,860 

-15.5 

-14.5 

-12.568 

-20.4 

-  19.0 

-21  1 

-22*6' 

OR  or 

-8.727 

-13,8 

-12.8 

EOPE 

-11,290 

-18,4 

-17.1 

-19.0 

“The  z  axis  is  the  principal  axis  of  the  molecule. 

'’Basis  set  1;  [5s3p2d/is2p]. 

''Basis  set  3;  [is3p2ii/3s2pld]  plus  two  s,p  lone  pair  functions  {exp.  0.1). 
'’Reference  50. 

^Reference  4. 


basis  set  3.  The  present  lone-pair  static  hyperpolarizabili¬ 
ties  are  close  to  the  static  values  obtained  by  Maroulis'^**^ 
with  a  large  (is6p5d2f /4s4p2d)  basis  set.  At  the  SCF  level 
his  ^  components  are  —9.4,  —7.6,  and  — 1.3  and  give  a 
of  —  1 1.0,  all  of  which  are  in  excellent  agreement  with  our 
results.  At  the  MBPT(2)  level  /3||  =  —  17.3,  compared  to 
our  —17.5  value.  In  a  somewhat  smaller  basis  {W\)  he 
obtains  a  fourth-order  triple  contribution  to  0||  of  —2.7 
a.u.  compared  to  our  CCSD(T)  measure  of  — 1.8  a.u.  It  is 
customary  for  fourth  order  to  overestimate  the  real  triple 
excitation  effects,  which  are  better  estimated  by  T(CCSD) 
or  the  (T)  correction.  The  theoretical  estimate  of  both 
hyperpolarizabilities  after  the  d’persion  correction  pro¬ 
vides  good  agreement  with  experiment  for  H2O. 

NH3 

Lazzeretti  and  Zanasi  (LN)’*  report  the  dipole  mo¬ 
ment,  polarizability  and  the  static  first  hyperpolarizability 
of  NH3  at  the  SCF  level  while  Liu  and  Dykstra  (LD)” 
also  report  SCF  static  first  and  second  hyperpolarizabili¬ 
ties.  The  LN  results  in  a  (4s6p3</l  f/6s2p\d)  basis  are  con¬ 
sistent  for  a  and  Aa  with  values  of  12.9  and  0.51 1,  respec¬ 
tively.  The  LD  results  are  12.64  and  0.55.  Maroulis 
reports  12.93  and  0.49  in  his  largest  basis 
(Is5p5d2f /Ss7>p\d)  which  are  in  excellent  agreement  with 
our  results  (Tables  XII  and  XIII).  Maroulis’  dipole  mo¬ 
ment  is  also  in  somewhat  better  agreement  with  ours,  being 
0.639  compared  to  our  0.637,  while  the  value  is  0.649  for 
LZ  and  0.642  for  LD.  We  see  a  more  dramatic  change  in 
)3i|  where  LN  get  12.86,  LD  14.84,  and  Maroulis  15.99, 
compared  to  our  SCF  basis  3  value  of  15.1.  For  the  indi¬ 
vidual  tensor  elements,  our  value  of  11.1  differs  sub¬ 
stantially  from  the  LN  value  of  7.8  and  the  LD  value  of 
9.1.  However,  it  is  close  to  the  M  1 1.8  value.  For  our 
basis  3  value  of  7.0  falls  among  the  LN  value  of  6.8,  the  LD 
value  of  7.8,  and  the  M  value  of  7.4. 


TABLE  XL  C'ofiipixients  ot  thv  p<>lari/-ahilit>  and  veumi  !iyj>fr(»'la/i/ 
abitities  of  ihc  H_.0  molecule  and  resuils  foi  \anous  pfatcccscc  im  a  u  ) 


SCF/ 

IDHF 

Mt)PI|2) 

CeSD 

cest)!  I  j 

1  .\P 

a„ 

7.8428 

7  9564 

9  6  .)62 

«>v 

9,1647 

10  0505 

9  8930 

9  1645 

10.0352 

9.8730 

10  0166 

Og, 

8  4993 

9  7503 

9  4881 

8  5264 

9  8193 

9.5416 

9  7285 

a 

8  50 

9.19 

9  55 

8  55 

9  87 

9  61 

9  19 

9  HI' 

\a 

1  14 

0.42 

0  56 

1  05 

0.25 

0  41 

0  34 

0  661' 

r,.,x 

1200.0 

2200  0 

20000 

1500  0 

2800  0 

27000 

29(8)  0 

470,0 

6300 

5900 

550.0 

820  0 

760.0 

8200 

y  rtrr 

770.0 

1200  0 

IIOOO 

9200 

1540  0 

1390  0 

1 5(X)  0 

Y.x,t 

330  0 

580,0 

550  0 

400.0 

7500 

700,0 

770.0 

2800 

4200 

400  0 

2800 

430,0 

4100 

4400 

290.0 

550.0 

540,0 

350.0 

660.0 

6200 

680  0 

Tii 

850.0 

1400.0 

13000 

lOlO.O 

1780.0 

1650.0 

1800  0 

THG 

1273.0 

21000 

2000.0 

1515.0 

2700  0 

25000 

2700,0 

dc-SHG 

1025.0 

1700.0 

16000 

1216.0 

2100,0 

2000,0 

2200  0 

2310  t  120* 

IDRI 

961.0 

1600.0 

1500.0 

1139.0 

2000.0 

1900  0 

2000.0 

EOKE  or 

904,0 

1500  0 

1400  0 

dcOR 

1069.0 

1900.0 

1700.0 

1900  0 

*M.  A.  Spackman.  J.  Phys.  Chem.  93,  7594  (1989),  as  obtained  from 
analysis  of  several  experiments. 

'’Reference  4. 


C.'.nsi.iering  y,  our  SCF  is  closely  bracketed  by 
that  of  LD,  1070  and  the  1297  value  of  M.  is  the  same 
to  two  digits  in  all  calculations.  Similarly  our  of  4500 
is  again  identical  to  two  significant  figures  with  M,  while 
being  only  slightly  larger  than  the  LD  value  of  4300. 

The  correlation  corrections  are  more  than  half  of  the 
total  hyperpolarizability,  emphasizinR  that  the  comparison 
to  experiment  at  the  SCF  level  is  meaningless.  The  only 
other  high-level  correlated  study  of  hyperpolarizabilities  is 
the  recent  work  of  Maroulis. '  In  comparisons  with 
his  SDQ(4)  results  in  a  (()s4p4d\f/4s2p\d)  basis,  he  ob¬ 
tains  =  30.47  which  is  nearly  identical  to  our  results  for 
CeSD,  the  infinite-order  extension  ofSDQ(4).  In  a  some¬ 
what  smaller  basis  he  evaluates  the  fourth-order  triples  to 
contribute  4.42  a.u.  to  P-  .  We  obtain  4.3  from  the  T  eval¬ 
uation  in  CCSD(T).  Hence,  all  our  results  for  pj  are  en¬ 
tirely  consistent  with  those  of  M.  Once  TDHF  dispersion  is 
added,  which  is  43%  for  SHG  (the  largest  we  have  ob¬ 
tained),  our  final  CCSD(T)  result  falls  within  the  experi¬ 
mental  error  bars.  The  other  correlated  result  is  limited  to 
second  order.'"  At  w  =  0,06  a.u.,  their  sectmd-order  SHG 
value  is  40  and  45.3  at  ira  =  0.07,  in  good  agreement  with 
the  present  results.  Similarly,  their  OR  value  is  .32.6  and 
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TABLE  XII  Dipiile  moment  and  components  of  the  first  hyperpiilan/ 
abilities  of  the  NH)  molecule  and  results  for  various  poKessts,  ( in  a  u  )  " 


SCF/IDHF 

MBPT(2) 

(  CSD 

CCSDl  T) 

F.XP 

0.6.164 

0.5994 

0.5974 

0.6369 

0.599! 

0  5976 

0  5901 

0.578‘ 

-6.4 

-  -10.2 

26  4 

-11.1 

-37.8 

-  33.'* 

.19,6 

P.;, 

-6.8 

-7.7 

-  7,6 

-7,0 

-  8.1 

-8.0 

8.8 

P 

-12  1 

-27.4 

-  25,0 

-15  1 

-32.5 

-  30.0 

-  .14,3 

SHG 

-  17.946 

-  41  0 

-37.3 

-21-978 

-46  6 

-42.5 

-49  1 

-  48  4a  1  2'’ 

OR  or 

-  13.595 

-31,0 

-28.3 

HOPE 

-  16,740 

-36.3 

-33.1 

-38.3 

*2  is  the  principal  axis  of  the  molecule. 

'’First  number  is  from  basis  set  1;  second  from  basis  set  3 
‘AV  H.  Weber,  J.  Mol.  Spectrosc.  107,  405  ( 1984) 

‘'Reference  4 

33,9  at  the  respective  frequencies.  As  the  comparable 
TDHF  percentage  dispersion  effect  is  virtually  identical  to 
that  obtained  in  a  second-order  correlated  dispersion  cal¬ 
culation,'*  our  final  CCSD(T)  results  should  be  an  excel¬ 
lent  estimate.  Vibrational  corrections  would  be  expected  to 
increase  the  magnitude  of  the  theoretical  values  somewhat. 
For  the  correlated  second  hyperpolarizability,  we 
again  have  consistency  between  the  present  basis  3  results 
and  those  of  Maroulis,  The  latter  SDQ(4)  result  of  3864  is 
again  nearly  identical  to  our  CCSD  38CX3  value.  Also,  the 
T4  estimate  in  a  smaller  basis  of  311  is  in  line  with  our 


correction  of  4{X),  so  that  Ins  esnrnaied  MBPl  fT  i  resuiis 
m  his  (bs4p4Jl//4.%2p]dj  basis  would  he  virtual!)  identi¬ 
cal  to  our  current  basts  3.  CCSD(  I  )  results  Including 
dispersion,  which  is  40'''/  for  dc-SH(J  results  we  obtain 
good  (89J-)  agreement  with  experiment  In  this  example, 
the  POL  basis  is  not  too  bad,  itself,  primarily  only  under¬ 
estimating  the  fi.,.  element  significantly  The  .VfHPTiJ) 
values  are  slightly  larger  than  CCSD  values  and  in  better 
agreement  with  experiment  for  both  the  HjO  and  NH, 
molecules.  It  should  be  noted,  however,  that  inclusion  of 
triple  excitations  in  general  brings  the  CCSD  hyperpolar- 
izabiiities  back  to  a  value  close  to  that  for  MBPT(2 ). 

HjS 

HiS  is  a  molecule  that  is  chemically  and  structurally 
quite  different  from  HvO,  It  involves  a  second-row  atom, 
and  has  an  HSH  angle  of  92.2°.  The  dipole  moment,  po¬ 
larizabilities,  and  hyperpoiarizabilities  of  HvS  are  summa¬ 
rized  in  Tables  XIV  and  XV  with  all  three  basis  sets.  The 
first  hyperpolarizability  is  one  of  the  smallest  among  the 
molecules  studied  here  and  therefore  a  very  careful  treat¬ 
ment  is  necessary  even  in  the  prediction  of  its  sign.  For 
example  the  SCF  values  obtained  here  have  an  opposite 
sign  to  that  of  experiment.  The  correlated  calculation  pre¬ 
dicts  the  same  sign  as  the  experiment,  but  rather  different 
values  for  different  basis  set  choices.  Once  again  the  dipole 
moment  and  polarizability  seem  to  be  adequate  in  any  of 
our  three  bases,  and  the  CCSD  is  slightly  belter  than 
that  given  by  CCSD(T).  However,  a  dramatic  change  is 
seen  in  the  hyperpoiarizabilities  in  the  results  between  the 


TABLE  XIII  Components  of  the  polarizability  and  second  hyperpoiarizabilities  of  the  NH,  molecule  and 
results  for  various  processes  (in  a.u  )  ’ 


SCF/aOHF 

MBPT{2) 

CCSD 

CCSD(T) 

EXP 

12  763 

13.779 

13.536 

12.7683 

13.7733 

13.5280 

13  7068 

a. 

13.279 

15,711 

15  242 

13,3475 

15  8334 

15  3346 

15  7139 

a 

12  94 

’4.42 

14.10 

12.96 

14.46 

14  13 

14  38 

14  56* 

Aa 

0.52 

1  93 

1,71 

0,58 

2,06 

181 

201 

1  94" 

1100.0 

1500.0 

14000 

1200.0 

ISOO.ri 

1700,0 

1800” 

Yaiz 

4200.0 

7300.0 

6900.0 

4500.0 

8000.0 

7400  0 

8200.0 

y XX2Z 

1100.0 

1800.0 

1700.0 

1100.0 

1900.0 

1800  0 

2000.0 

Til 

2300.0 

3700.0 

3500.0 

2400.0 

4100.0 

3800.0 

4200.0 

THG 

4590.0 

7400.0 

7000  0 

49250 

83000 

7800.0 

8500  0 

dc-SHG 

3072.0 

5000.0 

4700.0 

3276.0 

5500,0 

52000 

5600.0 

6090*  lift 

IDRI 

2762.0 

4500.0 

4200.0 

2937.0 

5000.0 

4600.0 

5100.0 

OKE  or  dcOR 

2301  0 

4100.0 

3800.0 

2652.0 

4500.0 

4200.0 

46000 

•G.  D,  Zeiss  and  W.  J  Meath,  Mol  Phys  33,  1155  (1977) 

’’N.  J  Bridge  and  A.  D.  Buckingham.  Proc.  R  Soc.  London  Ser.  A  295,  334  (1966). 
‘  Reference  4. 
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TABLE  XtV  Dipiilt;  nuimeiit  ai,v>  ^ompuncafs  of  (he  lirs!  livivrfKiian/- 
abililies  uT  the  tl;S  miileeule  and  results  for  sarious  pnxesses  im  a  u  )  ' 


SCF,  TDHF 

MBPT(2) 

CCSD 

CCSD ( 1  ) 

FXt> 

Me 

04.375 

0  4053 

0  1924 

0  4389 

0.4029 

0.39  j 

0  .1814 

0.4,184 

0  4037 

•  J 

0  3850 

0  401‘ 

U  1 

6  2 

\1 

6,8 

!.8 

0  8 

0  7 

7.9 

0.2 

2  6 

II 

A:.. 

8.2 

20 

■  12 

6.2 

5  5 

2.1 

-4.2 

-  5.5 

2  2 

4.4 

Av. 

6 

-  11.4 

1 1  7 

-  - 10,7 

10.2 

-  10  6 

10  6 

9,9 

-  9.2 

9.5 

,4  5 

0 

4  7 

-4  4 

-  1  7 

13 

-  10  5 

9  3 

2.2 

--8.8 

-  5.5 

-  7.7 

SHG 

5.813 

-5  5 

-  2  1 

1.488 

•  11.6 

--  7.9 

10.3 

2  380 

- 10.1 

-6.4 

-  8.8 

9  95  1-2.L 

OR  or 

5.038 

.,.4.8 

-  1  9 

EOPE 

1  450 

-  11,3 

-  7.7 

-  10  0 

2212 

-94 

-  5  9 

-  8  2 

‘Results  are  for  basis  1,  2.  and  3,  i  axis  is  principal  axis  of  the  molecule 
'’R  E,  Hilgcrand  M  W.  P  Strandberg,  Phys  Rev  83,  575  (  1951 1;  C  .A 
Bums,  }  Chem  Phys  28.  427  (1958),  G  A  Bums,  li/d  31.  127(3 
(1959), 

Reference  4 

POL  basis  1  and  basis  2,  due  to  the  diffuse  H  d  function  in 
the  latter  The  addition  of  the  lone  pair  functions  does  not 
necessarily  improve  the  answer  as  there  is  still  significant 
oscillation  between  basis  2  and  basis  3.  This  may  be  partly 
due  to  the  location  of  the  lone  pair  functions  which  is 
rather  different  than  for  HiO.  Also,  the  effect  of  the  addi¬ 
tional  electrons  and  the  comparative  degree  of  reliability  of 
the  starting  POL  basis  for  second-row  atoms  can  also  be  a 
consideration.  Basis  2  gives  better  agreement  with  experi¬ 
ment,  but  the  oscillatory  behavior  suggests  further  basis  set 
extensions  should  be  considered  before  claiming  conver¬ 
gence. 

Y  shows  a  more  modest  change  with  basis,  but  basis 
3  moves  the  result  in  the  wrong  direction  compared  to  the 
experimental  value.  There  seems  to  be  room  for  important 
basis  set  improvements,  although  generally  contrary  to  our 
experience  ( is  larger  with  every  basis  set  increase  in  this 
paper)  y  .  seems  to  be  overestimated  in  basis  3  compared 
to  the  dc-SHG  experiment. 

DISCUSSION 

In  Tables  XVI  and  XVII,  we  summarize  theoretical 
results  of  the  first-  and  second-electric  susceptibilities  for 
the  several  small  molecules  studied  in  this  paper,  now  in 
esu  units  as  reported  in  the  experimental  papers  (see  foot¬ 
notes  to  Tables  II  and  VI  for  conversion  factors).  The 
electric  susceptibilities  are  -  (5-, /2  and  =  y  /6. 
The  basis  sets  used  in  this  study  are  created  systematically 
from  polarizability-consistent  basis  sets^'’  and  are  only  of 
modest  size.  The  intent  is  to  provide  a  broad  study  for 
several  molecules  Hence,  no  particular  effort  is  made  to 
obtain  the  hyperpolarizahilities  of  a  specific  molecule. 


I  he  results,  taken  ax  a  whole,  allow  us  to  draw  some 
general  concliisions  about  hyperpoian/abiliiics  of  small 
molecules  For  y  ‘  the  static  SCF  resulis  have  an  average 
error  of  56^’f  compared  to  the  quoted  experirnciiial  values, 
s's  at  least  for  small  molecules  we  should  not  e.xpeci  St  I 
results  to  provide  better  than  a  factor  of  2  agreement  F di 
y“'  the  error  is  a  similar  43G  The  average  SHG  IDllf 
frequency  dependence  for  the  ruby  laser  frequency  of  2 1 '  f 
improves  v'''  results  to  a  4H'T  error;  i.e  .  for  SHG  the 
average  dispersion  is  about  of  experiment  The  average 
dispersion  for  y  ’’  dc-SHti  is  a  larger  2bG  of  ex¬ 

periment),  reducing  the  average  error  for  TDHI'  to  2H'^; 
Assuming  the  scaled  TDHF  dispersion  approximaiion 
used  throughout  this  paper,  the  addition  of  secenui-i  adc! 
correlation  reduces  the  y'"'  error  to  IIG  of  expcnrneiu, 
and  if  we  exclude  the  F'H  example,  whose  experimental 
value  we  question,  the  error  is  b%.  I'he  370^  effect  of 
correlation  on  y‘7*  is  about  tour  times  the  8'’)  dispersion 
effect.  For  y’’’,  the  .VIBPT(2)  corr.lc'ion  effect  is  i8G,  or 
slightly  larger  than  the  dc-SHG  dispersion  effeci,  reducing 
the  residual  error  for  y'"  to  lOG  Hence,  except  for  FH. 
in  a  reasonable  basis  TDHF  dispersion  weighted 
MBPT(2)  gives  an  average  error  of  lOG  for  y  '  The 
infinite-order  correlation  corrections  inircxluced  by  CCSD 
actually  diminish  the  effect  of  correlation  somewhat 
(  -8G)  for  y'’’  although  for  y'’'  CCSU  and  MBPT(2) 
have  the  same  average  error.  Once  triple  excitations  are 
included  via  CCSD(T),  the  average  error  for  TDHF 
scaled  dispersion  CCSD(T)  results  is  reduced  to  10. 7G  for 
y':''  and  10, IG  for  y'  '*.  Again  the  exclusion  of  FH  for  y  ‘ 
would  reduce  the  average  error  to  6,4G  for  y  '  Triple 
excitations  account  for  lOG  for  y  '  .  but  less  than  IG  on 
average  for  y*  Clearly,  even  larger  basis  sets  will  change 
the  distribution  of  effects  somewhat,  but  the  level  empk'y  ed 
in  the  present  calculations  should  be  sufficiently  high  to  be 
indicative  of  the  various  contributions.  Finally,  it  appears 
that  goixl  basis  set,  correlated,  frequency  dependent  results 
that  are  within  about  lOG  of  the  exfierimental  values  are 
possible.  Exclusion  of  any  one  of  the  above  contributions, 
however,  would  destroy  this  degree  of  agreement.  and 
f'.ih  values  are  quite  small  m  the  SHG  and  dc-SHG  cases.  '■* 
Assuming  adequate  basis  set  convergence,  the  leniatn- 
ing  primary  uncertainty  lies  in  the  inevitabie  limitations  in 
the  treatment  of  correlation.  For  the  current  results  ihis 
would  particularly  pertain  to  the  role  of  triple  excitations 
CCSD(T),  which  is  the  noniterative  approximaiion  to 
CCSDT-1,*^‘  and  has  been  found  to  mimic  well  full  CCSDT 
in  several  cases  for  electronic  energies  and  vibrational  fre¬ 
quencies.'^"’  As  such  it  IS  generally  better  than  the  original 
noniterative  CCSD  f  T( CCSD)  method,’''  which  is  equiv 
alent  through  fourth  order,  but  differs  by  a  single  fifth 
order  term,  particularly  in  ditTn  nil.  multircfcrcncc 

cases  None  of  the  current  molecules  arc  expected  to  ex¬ 
hibit  any  unusual  nuiltireference  character,  but  an  assess 
ment  of  reliability  of  CCSD{1  )  for  ckrlnc  propt-rlics  has 
not  yet  been  made.  Since  ITfCCSD);  is  norniallv  greater 
than  i(T)j,  when  the  two  triple  estimates  arc  close,  the 
extra  term  matters  little  and  we  usuallv  have  confidence  in 
the  triple  excitation  estimate  However,  for  some  c.iscs. 
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TABLE  XV  Cs'trifHinenl>  of  ihe  polanzabilny  aiid  sucutid  hypciTK'lari/abiliiifi  of  the  H_S  molecule  und 
results  for  curious  priK'essc's  (in  a,u, )  ' 


SCF/TDHK 

Mnt’r(2) 

CeSD 

C'CSOiT) 

KXC 

’T80SG 

25  1651 

24.8842 

2.V894! 

25  3435 

25  0502 

25  3644 

23  8.W8 

25.3537 

25  0353 

25.3507 

a 

23.6693 

24-2252 

23  99*)? 

24.3423 

24  1047 

24.2784 

2.3  7556 

24,4116 

24,1738 

24.3484 

(X., 

23.. >747 

24.234f) 

23.9832 

23-4207 

24.3808 

24  1187 

24  .3251 

2.3  430,3 

24.45.30 

24,1815 

24  3464 

a 

23  62 

24.54 

24  2^> 

23  67 

24.60 

24  43 

24  66 

23  68 

24  74 

24  46 

24  70 

25  5.''  26  6- 

Aa 

0  38 

0  44 

0  84 

U-41 

0  48 

0  44 

1  C)6 

0.37 

0.^2 

0.86 

0  44 

y  I  ix  t 

10  iy)o.o 

12900  0 

1 2  6(X).0 

10  lOftO 

13  000,0 

12  7(X).0 

1.3  5000 

10  900.0 

14  200.0 

!  3  500.0 

14  400  0 

2  100  0 

2  600  0 

2  500.0 

?  500,0 

3  200.0 

3  000  0 

^  0 

2  6(X)  ,0 

3  300.0 

3  21X1.0 

3  400.0 

3  2(X)  0 

4  000  0 

.3  800  0 

3  600.0 

4  7fXJ.O 

4  500  0 

4  7(X).0 

3  900.0 

5  100,0 

4  800.0 

5  KX)  0 

r,.,, 

2  .300  0 

3  (,XX).0 

2  900,0 

2  4('X),0 

3  3(X).0 

.3  100.0 

3  300,0 

2  6C00 

3  500.0 

3  .300.0 

.3  500,0 

V  .... 

1  3(XJ,0 

1  600  0 

!  600,0 

1  200,0 

1  500.0 

1  500  0 

1  500  0 

I  200.0 

1  600.0 

1  500.0 

1  600.0 

y<..> 

2  2fX),0 

2  800.0 

2  700.0 

2  400.0 

3  200.0 

3  1000 

.3  300  0 

2  400.0 

3  300.0 

3  too  o 

3  ,3000 

>' 

5  400  0 

6  900  0 

6  700  0 

5  700.0 

7  400.0 

7  1000 

7  500.0 

6  000.0 

7  8(X).0 

7  500.0 

7  9(X)0 

THG 

13  429  0 

17  200  0 

16  700.0 

14  290.0 

18  700,0 

18(X)0  0 

19  000  0 

15  223.0 

20  200.0 

19  3(X)  0 

20  000  0 

dc-SHG 

7  8150 

10  000  0 

9  700.0 

8  275  0 

10  800  0 

10  400  0 

1 1  exx)  0 

8713,0 

11  600,0 

1 1  oai.o 

1 1  700  0 

10  .300  a  260” 

IDRI 

6  840  0 

8  700,0 

8  500.0 

7  221.0 

9  500.0 

9  100  0 

9  600  0 

7  573,0 

10  1000 

9  6(X)  0 

10  200.0 

F:OKEor 

6  037.0 

7  700.0 

7  500  0 

dcOR 

6  360.0 

8  300,0 

8  (XX).0 

8  5000 

6  650.0 

8  800,0 

8  400.0 

8  900  0 

■‘Results  are  for  basis  1.  2,  and  3. 

’’J  O  Hirschfeldcr.  C.  F.  Curtis,  and  R.  B.  Bird.  Molecular  Theory  of  Casef  and  Liquids  ( Wiley,  New  York, 
1054).  p.  950. 

A  A  Maryotl  and  F  Buckley,  Nat,  Bur.  Stand.  (U  S  )  Circ.  537  (  195.3) 

’Reference  4, 


though  CCSD  +  T(CCSD)  and  CCSD(T)  coincide,  both 
can  differ  significantly  from  CCSDT."*^  In  Table  XVIII  we 
compare  CCSD  f  T(CCSD)  and  CCSD(T)  for  static  hy¬ 
perpolarizabilities. 

For  all  cases  but  FH,  the  y.i  values,  retaining  two 
significant  digits,  are  essentially  the  same.  For  the  av¬ 
erage  increase  in  magnitude  is  about  7%  except  for  CO, 
which  shows  almost  no  change.  Hence,  we  are  confident  of 
the  triple  estimate  for  these  examples.  Some  small  change 


would  occur  in  comparisons  of  /i.  with  experiment  for 
FH,  HyO,  and  NH,  if  we  used  T(CCSD),  but  since  we 
tend  to  believe  that  T(CCSD)  overestimates  the  triple  cor¬ 
rection,  we  prefer  the  CCSD(T)  values.  However,  in  the 
absence  of  full  Cl,  or  results  from  higher  approximations 
like  CeSDT*^  or  CCSD'i  Q,'’'  we  can  be  relatively  (but  not 
absolutely)  certain  that  the  triple  excitation  estimate  in 
CCSD(T)  will  be  reliable  for  quantities  as  sensitive  as  hy 
perpolarizabilities.  Here  except  for  the  very  small  liyper- 
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TABLE  XVI  Thwretical*  and  expenmenlal*'  hyperpulanzabiliiies 
( 10  ''  esu/molecule)  for  SHO 


MBPT(2) 

,yP’ 

CCSD 

CCSD(T) 

Exp 

CO 

II  2 

11.4 

11.7 

12  9*1  4 

HF 

-3,3 

-3.2 

-3.4 

-4  70  *  0.41 

H,0 

-8.8 

-8,2 

-9  1 

-94*04 

NH, 

-20.1 

-18.4 

-21.2 

-20  9  *  05 

HjS 

-4.4 

-2.8 

-3.8 

-4,3*09 

'Value  corrected  for  the  dispersion  effect  at  694.3  nm  using  the  TDHF 
percentage  dispersion  correction.  All  the  calculations  are  performed  with 
basis  sets  (5s.V2£f]  for  C.  N,  O,  and  F;  I7s5p2d]  for  S;  and  [Jslpld]  for 
H.  The  lone-pair  functions  are  added  for  HF,  H,0,  NHi,  and  H^S  All 
molecules  at  experimental  geometries  and  there  is  no  estimate  of  vibra¬ 
tional  corrections. 

'Value  obtained  by  dc-induced  second  harmonic  generation  (Refs.  4,  5, 
and  56). 

polarizabilities  of  FH,  any  difference  with  experiment  due 
to  triples  estimates  will  not  affect  our  10%  error. 

Because  of  cancellation  of  errors,  MBPT(2)  can  pro¬ 
vide  results  for  correlation  that  can  be  in  accidentally  bet¬ 
ter  agreement  with  experiment  than  the  theoretically  more 
complete  CCSD  as  the  effect  of  triple  excitations  is  gener¬ 
ally  to  move  CCSD  closer  to  MBPT(2).  This  computa¬ 
tionally  useful  observation  is  not  universal,  nowever,  as  the 
C2H4  example  shows  a  significant  MBPT(2)  overestimate 
for  compared  to  CCSD(T). 

Because  of  the  above  conclusion  that  calculations  at 
the  level  used  here  provide  results  to  within  a  mean  error  of 
10%,  the  FH  example  still  appears  anomalous,  as  it  did  in 
our  earlier  work.^  ®  In  basis  3,  at  the  CCSD(T)  error  is 
28%  for  and  24%  for  Although  not  too  great, 
they  are  about  three  and  two  times  the  next  largest  percent 
error,  respectively.  Compared  just  to  H2O  and  NHy,  the 
FH  error  is  close  to  a  factor  of  10  for  and  3  for 
Of  course,  FH  can  also  have  characteristics  that  are  pecu¬ 
liar  to  it,  which  selectively  reflect  deficiencies  in  the  ap¬ 
proximate  quantum  chemical  solutions  employed,  and  any 
systematic  error  would  be  amplified  in  the  percentage  cor- 


TABLE  XVII.  Theoretical'  versus  experimental’’  hyperpolanzabilities  of 
the  molecules  { 10  *  ”  esu/molecule). 


MBPT(2) 

t'” 

CCSD 

CCSDtT) 

Exp- 

58.2 

59.3 

59  3 

60  5  (exact) 

CO 

151,0 

149.0 

160.0 

144  0  *  4,0 

N, 

88.7 

90.6 

96,0 

86  6*  1.0 

CO2 

Ills 

107,9 

1 100 

1119*  1.3 

C,H. 

960  0 

820,0 

860.0 

758.0*  17  0 

HF 

520 

49.0 

53.0 

70.0*  100 

HjO 

180  0 

170.0 

180.0 

194.0*  10.0 

NH, 

460.0 

4.300 

470,0 

5110  *  9.0 

,>  i 

9700 

9.300 

980.0 

865.0*22  0 

'All  the  calculations  are  performed  with  basis  sets  [5r}p2d]  for  C,  N.  O. 
and  F;  (7,v5p2</)  for  S;  and  (3j2pld)  for  H  The  lone-pair  functions  are 
added  for  HF.  HjO,  NH,.  and  H^S  All  molecules  at  experimental  ge¬ 
ometries  and  there  is  no  estimate  of  vibrational  corrections 
’’Value  obtained  by  dc-induced  second  harmonic  generation  (Refs  4.  5, 
and  56) 


TABLE  XVlll  Comparison  of  diH  crciiT  nnnilcratt'»c'  fnplc 
estimates  for  static  hyi'nrrpi’ilan/abfiiiies  ut  u  > 


CCSDiT) 

n 

C  CSD  -  If  CCSD) 

CCSDi  1  ! 

CC.'S1>  •  KCCSD) 

CO; 

lUX) 

i 

N. 

IfMX) 

HUM) 

C,H, 

67(X) 

67(X) 

CO 

23  5 

230 

16tX) 

IWX) 

HF 

-7  3 

-  7  8 

5<)0 

ft!U 

H,0 

-  18  0 

-  19,5 

1 8(X) 

!9(Xj 

NH. 

-34  3 

3ft  9 

42fX) 

43(.K) 

H;S 

-  7  7 

8  3 

79(X) 

rection  since  it  has  the  smallest  v  '  ■  ‘’f 

three  molecules.  In  other  words,  systematic  errors  that  af¬ 
fect  all  calculations,  will  disproportionately  cause  a  higher 
p»ercentage  error  in  a  small  value  than  in  a  large  one. 
Hence,  the  molecule  requires  further  consideration. 

The  original  paper  on  FH  (Ref.  6)  used  the  SDQ- 
MBPT(4)  method  (the  fourth-order  approximation  to 
CCSD)  and  computed  static  hyperpolanzabilities  in  sev¬ 
eral  different,  large  polarized  basis  sets.  The  D  {i>s5p4(I/ 
5s2p)  basis  gave  —2.83  .<10  ”  esu/molecule  for  v  * 
((i)  =  0)  at  and  39x  10  esu/molecule  for  y’  ’’  (f/j  =  0). 
These  numbers  should  be  compared  with  the  current  basis 
3,  CCSD  static  values  which  are  —2,9  and  44,  respectively 
In  other  words,  even  the  rather  drastic  and  unconventional 
ba.sis  set  addition  of  lone  pair  functions  has  only  had  a 
modest  effect  on  the  original  1979  values.  Comparison  with 
the  later  1986  results,®  namely,  —2.79  and  37,  in  basis  E 
{i.e.,  (6s5p4d2  //533/7)]  are  again  completely  consistent,  In 
basis  F.  we  added  a  set  of  d  functions  on  H.  exponent  1 .0, 
to  the  E  basis  to  gne  {6s5pdd2  f/5sip\d).  Since  we  did 
not  previously  report  correlated  calculations  in  that  basis, 
we  do  so  here.  At  the  CCSD  level,  basis  F  gives  -  2.8  and 
37  esu  for  10*’  and  lO'®  y'  ”,  respectively,  and  they  are 
again  consistent.  The  F  bass  MBPT(2)  values  are  -  2.9 
and  41.  Further  mi,dificatiOii  of  bases  E  and  F  led  to  our 
basis  I,®  namely.  (6s5p5oi  j /5s?ip),  which  had  more  diffuse 
d  and  /  exponents.  This  basis  gave  a  larger  /3«z  clement, 
and  the  functions  in  the  basis  have  since  been  found  to  be 
important  in  the  satisfaction  of  the  dipole  length  and  ve¬ 
locity  equivalence  in  RPA  calculations.'^  We  previously 
only  reported  MBPT(2;  calculations  with  basis  !.  Here, 
we  take  basis  I  and  add  the  d  function  on  H  ( exponent  1 ,0) 
to  give  what  we  will  call  basis  J  for  consistency  with  the 
prior  calculations,  and  do  a  full  set  of  correlated  calcula¬ 
tions.  In  basis  J  CCSD  values  are  2.97  and  40.  and  the 
MBPT(2)  values  3,08  and  43.  The  conclusion  is  that  de¬ 
spite  a  variety  of  conventional  and  unconventional  bases, 
our  current  static  results  differ  rather  insignificantly 
(  <  10% )  from  the  original,  static  r ,  values. 

The  first  additional  effect  that  more  significantly 
changes  the  computed  static  values  is  the  inclusion  of  triple 
excitation  effects.  This  contribution  in  basis  E  causes  the 
SDTQ-VIBPT(4)  values  to  be  3.13  and  43.  respectively. 
This  is  similarly  observed  in  the  current  basis  3  results, 
whose  values  at  the<’CSD(T)  level  are  3  15  and  47.  The 
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TABLE  XIX  Comparison  of  ekclric  susceptibility  results  for  I  H  from  difTcrcnt  calculations  (csu/molcculci 


Reference 

Level 

Basis 

lO'-'v'-'tw  - 0) 

10'%  ‘  (w  0) 

lO  'i  ' 

( 2, .M,, 

1979  BP  (6) 

SDg-MBPT(4) 

£);(6v5p4<//5s.V) 

-  2  83 

W' 

’  3  ( 

.3.5  r 

44!  44  ) 

.V(M 

1  3  20)“' 

44i45i'' 

MBPT(2) 

-2.73 

40” 

1986  SB  (8) 

CCSD 

f  (5s5p4rf2  //5.-.V) 

-  2  62 

35 

CCSDrT4 

2  97 

41 

.1  14 

(  3  29,)'' 

-14(4Si' 

MBPT(2) 

-2.72 

SDQ-MBPT(4) 

-2.79 

M 

MBPT(4) 

" 

3  13 

43 

3  30 

(  3  45  )' 

46(47  )•’ 

MBPT(2) 

/:(6j5p5r/l//5.v3/)) 

-  3.2 

41 

(CCSD  I-T4) 

-3  5 

44 

-  3,6  ( 

.  3  8)' 

47(48)' 

1990  CL  (14) 

MCPF 

(8.s«;>4J3//5i3/)2z/) 

-3.12 

43 

CCSD 

-2.89 

CCSD(T) 

-  3. 10 

43 

1992  JJJ  (28) 

MCLR-CAS  4220” 

/isipld)* 

-2  80 

37 

-3.08 

(  .  3 Tb)' 

4.3 

MCLR-CAS  6331 

-2.87 

1992  SB  (present  work) 

CCSDCT) 

J  (6s5p5d[  f/S’ipld) 

-3.21 

44 

.  .  3  49 

(-  3  65)' 

49(50)“’ 

CCSD 

" 

-2.97 

40 

MBPT(2) 

-3  07 

43 

CCSD(T) 

{isipld/is2p\d 
//3x(1.5,lp)  lone 

-3.15 

47 

-  3.43 

(  -3.59) 

5,3(55) 

pair) 

MBPT(2) 

-3.01 

47 

CCSD 

-2.93 

44 

CCSD(T) 

Largest  effects* 

-3  26 

47 

-3  58 

(-3.86) 

56(57) 

Experimenl 

(  -4.70 ±0.41) 

(70±  10) 

*(u  =  0.0656  a.u. 

'’Reference  56  reports  35  by  estimating  Footnote  38(a)  in  the  reference  corrects  this  value  to  39  using  our  computed  value  F{  !r  )  core  correlation 
effects  included. 


'First  line  incorporates  14%  dispersion  and  10%  zero-point  correction  as  estimated  by  Dudley  and  Ward  ( Ref  56),  while  the  second  value  was  obtained 
using  actual  TDHF  percentage  dispersion  and  computed  vibrational  average. 

‘'Vibrationally  "averaged"  results  in  {  )  Reference  31  finds  a  -0.28  unit  effect  for  y'’'(  -  Imo.o)  due  to  zero-point  average.  SB  found  -  0  15  at  the 
static  CCSD-(-T4  level  in  Ref  8  and  1  unit  for  y';'’(0',0.0.0).  In  the  new  results  we  report  values  at  ri,=  1  7621  a  u  (Ref  50)  instead  of  r,  to  partly 
incorporate  this  effect,  getting  -2.52  SCF,  -  3. 14  MBPT(2),  —  3.09  CCSD.  and  —3.31  for  CCSD(T).  Similarly,  the  r„10'''y'  '  (0:0.0,01  values  are  28 
SCF,  48  MBPT(2),  45  CCSD,  and  49  for  CCSD(T).  These  results  are  nearly  identical  to  our  previous  zero-point  averaged  corrections  from  l^Sb  ( Ref 
8). 

This  CCSD-cT4  result  from  Ref  8  was  extrapolated  from  the  MBPT(2)  results  for  basis  I.  assuming  all  effects  beyond  MBPT(2)  to  be  additive 
^This  ANO  basis  (Ref  31)  differs  from  that  of  Chong  and  Langhoff  (Ref  14). 

^We  report  results  from  the  largest,  computed  static  tensor  elements  from  any  of  the  above  calculations  (sec  the  text)  augmented  with  the  largest 
observed  zer(,-(X)int  vibration  and  dispersion  effects  from  Jaszunski  ei  al  (Ref  31 ).  The  MCLR  SHG  dispersion  is  10%,  and  di-SHG  dispersion,  19% 
This  would  seem  to  offer  an  upper  bound  to  the  magnitude  of  and  y'"  as  can  be  justified  by  alt  current  high-level  correlated  calculations 


next  important  element  is  to  average  over  the  zero-point 
vibration.’  Here  we  approximate  this  average  by  evaluation 
at  rj).  With  displacement  becomes  — 10.3  a.u.  com¬ 
pared  to  —9.6  a.u.  at  while  is  —  1.26  a.u.  instead  of 
-1.27  a.u.  at  r^.  For  CCSD(T)  this  changes  iO’’  y||”  to 
—  3  31  and  10^'’  to  49.  Finally,  using  the  TDHF  per¬ 
centage  dispersion,  we  obtain  -3.59  and  55,  respectively. 
If  we  go  back  and  compare  with  our  original  1979  results, 
including  the  upper  bound  10%  zero-point  vibration  cor¬ 
rection  and  14%  dispersion  correction  estimated  by  Dud¬ 
ley  and  Ward,^^  we  have  —3.5  and  49,  which  is  still  in 
good  agreement.  Alternatively,  we  can  use  the  actual  vi¬ 
bration  and  dispersion  correction  we  obtained**  to  give  the 
second  set  of  values  in  Table  XIX,  which  are  slightly 
smaller  in  magnitude.  So  despite  improvements  in  the  cal¬ 
culation  primarily  due  to  the  addition  of  triple  excitations, 
the  original  values  are  still  close  to  the  best  that  we  have 
since  obtained.  Hence,  our  current  values  seem  to  be  well 
converged. 

This  conclusion  is  also  supported  by  the  other  high- 
level  correlated  calculitions  that  have  been  reported. 


Chong  and  Langhoff***  use  a  more  extensive  atomic  natural 
orbital  (ANO)  basis'*^  (see  Table  XIX)  and  present  a 
modified  coupled  pair  functional  (VfCPF)  calculation  in 
addition  to  CC  results.  Their  static  values  of  —3.12  for 
lO'^  43  for  10*“  are  in  excellent  agreement 

with  all  prior  results  that  include  triple  excitations.  Also, 
the  CCSD(T)  results  are  in  excellent  agreement  with  the 
MCPF  values.  As  the  ANO  basis  is  generated  from  a  dif¬ 
ferent  premise  than  the  basis  sets  introduced  in  our  work, 
the  correspondence  of  results  is  indicative  of  a  large  degree 
of  basis  set  convergence.  Even  the  potentially  drastic 
changes  that  are  introduced  via  lone  pai-  functions  sh  'w 
little  difference  from  the  other  results,  contributing  a  max¬ 
imum  change  of  about  0,1  esu/molecule  in  10*’  v’’’  and  4 
esu/molecule  in  lO'**  y''*. 

The  other  high-level  calculation  reported  for  FH  is 
that  of  Jaszunski  eta/,'*  These  authors  report  MCLR- 
CASSCF  results,  Their  static  correlated  values  are  smaller 
in  magnitude  than  those  obtained  by  CCSD(T)  methods, 
possibly  due  to  the  neglect  of  residual  dynamic  (i.e  .  CD 
correlation.  Their  CASSCF  results  use  very  large  CAS 
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configuration  spaces  of  up  to  125  000  determinants  and 
perturbed  CASSCF  orbitals.  In  comparison,  our  J-basis 
CCSD(T)  explicitly  involves  singles  and  doubles 
( 144  000  determinants  and  implicitly  most  of  the  effect 
of  triples  and  quadruples!)  with  relaxed  (i.e.,  coupled  per¬ 
turbed)  Hartree-Fock  orbitals.  Their  CASSCF  results  are 
close  to  our  previous  static  SDQ-MBPT(4)  and  CCSD 
values  and  the  CCSD  values  of  Chong  and  Langhoff,'"*  but 
triple  excitations  enhance  the  prior  values  somewhat.  Add¬ 
ing  their  zero-point  vibrational  average  and  dispersion, 
their  final  results  remain  quite  close  to  those  we  obtained  in 
basis  £■;*  but  are  somewhat  smaller  in  magnitude  than  the 
best  current  values.  It  should  be  noted  that  even  very  large 
CAS  spaces  do  not  necessarily  introduce  all  appropriate 
dynamic  correlation  effects,  frequently  requiring  the  inclu¬ 
sion  of  Cl  on  top  of  their  CASSCF.  This  is  reflected  in 
their  computed  polarizabilities  which  have  only  60%  of  the 
correlation  effect.  Any  such  dynamic  correlation  error  in 
j3i|  oryji  (and  one  would  expect  more  in  the  latter  because 
like  a  it  increases  with  molecular  size),  would  modify  their 
results  to  be  close  to  the  magnitude  of  the  current,  largest 
values.  As  it  is,  they  obtain  nearly  the  basis  set  E 
CCSD-)-T4  results  of  our  paper.*  In  that  paper,  we  state 
that  we  believe  the  magnitude  of  will  fall  between  our 
computed  value  of  3.3  and  our  maximum  estimated  value 
of  3.8.  Their  result  of  3.36  falls  within  the  specified  range. 
It  is  likely  that  additional  dynamic  correlation  will  increase 
the  magnitude  of  their  computed  values  to  be  more  in  line 
with  the  current  CCSD(T)  results  of  3.6  and  55. 

To  alleviate  suspicion  regarding  any  anomalous  behav¬ 
ior  of  the  lone  pair  basis  set,  we  also  report  new  CCSD(T) 
atom-centered  basis  results  for  basis  set  J.  Previously,  we 
recognized  that  the  tensor  component  for  FH  had  been 
found  to  be  larger  in  some  of  our  other  basis  calculations 
than  in  our  E  basis  results.*  This  raised  the  question  of  a 
potential  basis  set  deficiency  in  the  perpendicular  direction. 
To  investigate  this  we  considered  several  basis  modifica¬ 
tions  that  led  to  the  I  basis  estimate  made  in  our  earlier 
paper*  whose  MBPT(2)  value  for  was  —  1.6  compared 
to  —0.7  in  basis  E.  As  the  component  remained  nearly 
the  same,  assuming  the  additivity  of  other  effects,  the 
larger  MBPT(2)  value  for  P^^  in  basis  I  led  to  a  final, 
maximum  estimate  of  —3.8  esu/molecule  for  10^^  To 
better  document  the  effects  of  this  basis,  we  now  report 
correlated  results  in  basis  J,  which  is  I  plus  the  set  of 
Cartesian  d  (1.0  exponent)  functions  on  H,  i.e., 
(6j5p5(/l//553pl<f).  Though  this  is  not  that  extensive  a 
basis,  our  static  CCSD(T)  result  of  —3.21  is  close  to  the 
lone  pair  basis  results,  while  the  value  of  44  is  slightly 
lower.  After  dispersion  and  vibration  correction  the  results 
are  —3  65  and  50  for  lO’^  and  lO”  esu/molecule. 
At  a  given  level  of  correlation,  usually  the  more  complete 
the  basis  the  larger  will  be  Yw  ■  If  this  applies  here,  the  lone 
pair  basis  might  be  considered  to  be  superior. 

It  is  informative  to  estimate  an  upper  bound  for  the 
magnitude  of  and  For  Chong  and  Langhoff 
get  —1.38  a.u.  at  the  MCPF  level  while  Jaszunski  et  al. 
find  its  value  to  be  --  1.10  a.u.  with  CAS  6331,  compared  to 
our  current  basis  3  value  of  —1.27  a.u.  and  basis  J 


CCSD(T)  value  of  -  1,47  a  u.  The  other  component,  /i 
is  -9.29  a.u  for  MCPF,  K.H5  for  the  CAS  6.Ctl.  and 
—  9.43  a.u.  for  basts  J.  Our  current  basts  3  CCSD(  1)  /I.., 
of  —9.62  a.u.  has  the  largest  magnitude  Selecting  all  these 
upper  limit,  correlated  values  for  the  components  of  fi  and 
y,  the  largest  latter  components  are  our  current  basts  .1 
values  shown  in  Table  IX  [for  comparison  the  J  basts  val¬ 
ues  are  y,„,=  383,  6(X).  161.  2(X)  (2!0 

in  basis  3)]  together  with  the  largest  vibrational  average 
and  dispersion  correction  taken  from  the  MCLR  results, 
we  obtain  the  final,  estimated  CCSD(T)  result  tn  Table 
XIX.  This  procedure  would  seem  to  offer  an  upper  bound 
to  the  magnitudes  of  y'/'  and  y'  "  that  can  be  justified  by 
all  calculations.  Obviously,  such  an  estimate  should  not  be 
considered  to  be  superior  to  calculations  that  have  in¬ 
cluded  all  effects  consistently,  as  there  are  mutual  interac¬ 
tions  within  a  given  basis  and  level  of  calculation  that  will 
result  in  a  balance  in  the  sizes  of  tensor  components  that 
should  certainly  be  more  appropriate  in  a  rigorous  calcu¬ 
lation  than  taking  the  largest  individual  components.  Sim¬ 
ilarly,  our  smaller  correlated  vibrational  correction  is  not 
necessarily  less  accurate  than  the  MCLR  value  as  the  ab¬ 
sence  of  Cl  can  be  important:  nor  is  the  larger  MCLR 
dispersion,  though  more  justified,  necessarily  the  most  re¬ 
liable.  Nonetheless,  we  obtain  final  limit  values  of  —  3.9 
and  57  for  the  electric  susceptibilities.  This  estimated 
bound  still  shows  about  a  17%- 19%  error  compared  to  the 
center  of  the  experimental  values.  However,  taking  the 
lower  limits  from  the  error  bars,  we  would  have  -  3.9 
compared  to  -4,3  and  57  compared  to  60,  which  would 
fall  within  the  10%  error  we  would  expect.  Obviously, 
however,  this  seems  to  be  the  maximum  magnitude  that 
can  possibly  be  justified  by  any  calculation  to  date.  All 
“consistent"  calculations  give  smaller  magnitudes  leaving 
unresolved  the  —20%  error  for  FH.  In  particular  our  best 
calculation  gives  —3.6  and  55  in  the  lone  pair  basis  Our 
best  atom  centered  basis  result  is  that  in  basis  J,  resulting 
in  —3.6  and  50.  Much  larger  basis  calculations  are  reason¬ 
able  and  should  be  made;  but  in  any  event,  we  conclude 
that  theory  supports  a  value  of  — 3.6±0.3x  l0  ''  esu/ 
molecule  for  yi,^*  and  m  y'  - ’,  55  ±  5  X  10  esu/molecule. 
Finally,  if  we  a.ssume  that  CCSD(T)  does  not  reliably 
measure  the  correct  value  of  triples,  and  in  particular,  that 
the  usually  less  reliable  T(CCSD)  is  superior  in  this  case, 
we  would  add  —0.16  esu  units  to  y':‘'  and  3  esu  units  to 
yi|”,  still  within  the  specified  error  bars. 

CONCLUSIONS 

SCF  theory  predicts  the  magnitude  of  the  hyperpolar¬ 
izabilities  to  be  smaller  in  comparison  with  SHG  and  dc- 
SHG  experimental  values  for  all  molecules  studied  here 
typically  by  a  factor  of  2,  Both  the  dispersion  correction  at 
the  Hartree-Fock  level  (TDHF)  and  the  correlation  cor¬ 
rection  [MBPT(2),  CCSD,  or  CCSD(T)]  at  the  si  'ic 
level  increase  he  magnitude.  The  theoretical  values  at  the 
MBPT(2)  correlated  level  augmented  by  the  percentage 
dispersion  correction  evaluated  by  TDHF  show  good 
agreement  with  the  dc-SHG  experimental  values  for  most 
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of  the  molecules.  Overall  agreement  of  the  theoretical  hy¬ 
perpolarizabilities  with  experiment  is  encouraging  for  the 
application  of  this  level  of  theory  to  larger  systems.  The 
polar  molecules  which  have  lone  pairs  and  hydrogens  are 
likely  to  be  more  difficult  to  describe  well,  as  the  basis  set 
requirement  seems  to  be  more  severe.  At  a  minimum,  at 
least  one  set  of  diffuse  polarization  basis  functions  must  be 
added  to  the  original  polarizability  consistent  basis  set  for 
hydrogen,  while  additional  functions  must  be  added  to  the 
heavy  atom  basis.  Because  of  H,  vibrational  effects  are  also 
likely  to  be  more  important  for  these  molecules,  and  results 
should  be  corrected  for  that  effect.  By  including  the  lone- 
pair  functions,  the  comparatively  high  polarization  of  the 
lone-pair  electrons  can  be  described  more  properly.  Alter¬ 
natively,  a  large  number  of  diffuse  polarization  functions 
including  /.  g,  and  h  ■  orbitals  located  on  the  heavy  atom 
would  be  expected  to  perform  the  same  task.  The  atom 
centered  basis  would  offer  more  opportunity  to  extrapolate 
to  the  basis  set  limit  than  the  lone  pair  basis  procedure. 
The  latter,  however,  has  the  effect  of  potentially  eliminat¬ 
ing  any  residual  prejudice  in  standard  basis  generation 
methods.  Consistencies  between  results  in  both  types  of 
bases  recommends  good  (but  not  definitive)  basis  set  con¬ 
vergence.  At  a  given  level  of  theory,  all  basis  set  improve¬ 
ments  lead  to  larger  y.  values. 

Correlation  effects  are  important  for  all  the  molecules 
studied.  MBPT(2)  tends  to  provide  slightly  larger  values 
for  both  hyperpolarizabilities  and  fortuitously  better  agree¬ 
ment  with  experiment  than  does  CCSD  itself  However, 
the  triple  excitation  effect  as  measured  by  CCSD(T)  sig¬ 
nificantly  improves  upon  the  CCSD  predictions  for  the  0 
hyperpolarizability,  with  results  usually  moving  somewhat 
closer  to  those  obtained  by  MBPT(2).  This  suggests  that 
MBPT(2)  may  offer  an  adequate  measure  of  correlation 
for  large  molecule  applications.  Note,  however,  that  for 
ethylene,  the  MBPT(2)  estimate  is  too  large.  This  may 
indicate  that  MBPT(2)  might  be  insufficient  for  larger  de¬ 
localized  systems,  such  as  linear  polyenes. 

Since  any  potential  systematic  error  in  a  calculation  is 
magnified  when  the  experimental  value  is  small,  the  small 
values  for  the  hyperpolarizabilities  of  the  FH  molecule 
combined  with  the  potential  for  a  systematic  error  make  it 
the  most  difficult  case  among  the  molecules  studied,  ft 
seems  that  a  discrepancy  of  about  20%  between  theory  and 
experiment  remains.  Unless  some  unusual  feature  encoun¬ 
tered  in  further  investigation  of  basis  sets,  correlation,  dis¬ 
persion  corrections,  or  other  elements  emerge,  it  seems  un¬ 
likely  that  the  remaining  —20%  error — as  small  as  it  is — 
can  be  resolved.  Experimentally,  any  population  of  excited 
vibrational  state.-.'^'^^  in  the  experiment  would  dramatically 
increase  |  and  jyij  |  and  could  account  for  the  differ¬ 
ence.  Another  experimental  measurement  would  be  wel¬ 
come. 
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